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Abstract 

We review different aspects of field theory at zero and finite temper- 
ature, related to the theory of phase transitions. We discuss different 
renormalization conditions for the effective potential at zero tempera- 
ture, emphasizing in particular the MS renormalization scheme. Finite 
temperature field theory is discussed in the real and imaginary time 
formalisms, showing their equivalence in simple examples. Bubble nu- 
cleation by thermal tunneling, and the subsequent development of the 
phase transition is described in some detail. Some attention is also de- 
voted to the breakdown of the perturbative expansion and the infrared 
problem in the finite temperature field theory. Finally the application 
to baryogenesis at the electroweak phase transition is done in the Stan- 
dard Model and in the Minimal Supersymmetric Standard Model. In all 
cases we have translated the condition of not washing out any previously 
generated baryon asymmetry by upper bounds on the Higgs mass. 
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1 The effective potential at zero temperature 

The effective potential for quantumJield theories was originahy introduced 
by Euler, Heisenberg and SchwingerEJ, and appHed to studies of spontaneous 
symmetry breaking by Goldstone, Salam, Weinberg and Jona-LasinioH. Calcu- 
lations of the effective_potential were initially performed at one-loop by Cole- 
man and E. ^einbergB and at higher-loop by JackiwQ and Iliopoulos, Itzykson 
and Martin 0. More recently the effective potential has been the subject of a 
vivid investigation, especially related to its invariance under the renormaliza- 
tion group. I will try to review, in this section, the main ideas and update the 
latest developments on the effective potential. 



1.1 Generating functionals 

To fix the ideas, let us consider the theory described by a scalar field (p with a 
lagrangian density C{(j){x)} and an action 

= J d^xC{cl,{x)} (f) 

The generating functional (vacuum-to- vacuum amplitude) is given by the path- 
integral representation. 



Z[j] = (Oout I Oi„)j = J d(^exp{z(5[(/.] + <j>j)} (2) 

where we are using the notation 

0j = J d^xc^{x)j{x) (3) 

Using (||) one can obtain the connected generating functional W[j] defined as, 

Z[j\ = exp{zM/[j]} (4) 
and the effective action r[0] as the Legendre transform of (^) 

m = wm- j d''x^-^j{x) (5) 

where 

■7,( ^ ^^^^ 
bj[x) 
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In particular, from and (^), the following equality can be easily proven, 



<5r[0] _ 6W[j] Sg _ . 



(7) 



5<l) 5j 

where we have made use of the notation (||) . Eq. implies in particular that 



(8) 



3=0 



which defines de vacuum of the theory in the absence of external sources. 

We can now expand Z[j] {W[j]) in a power series of j, to obtain its repre- 
sentation in terms of Green hmctions (connected Green functions G(^nj) 
as, 



(9) 



and 



Similarly the effective action can be expanded in powers of (fi as 

- °° 1 /■ - - 

- / d^xi...d^x^cj){x,)...c^ixn)r'^^Hx,,...,x^) (11) 



where F^") are the one-particle irreducible (IPI) Green functions. 
We can Fourier transform F*^") and (j) as. 



/7l 
n 
i=l 



d^ 



exp{ipja;i} 



Up) ^ I d^xe-'P'=(j>{x) 
and obtain for ^ the expression, 

d'^Pi 2 



i2n f S^^Hpi + ■■■+ p,,)T'-^\p) (12) 

(13) 



oo ^ n 

E/n 



(27r) 



iH-Pi) 



{2ti)H^^\pi + • • • +p„)F(")(pi, . . . ,p„) (14) 



n=0 i=l 

In a translationally invariant theory, 

(j){x) = (t)c 



(15) 



the field (j) is constant. Removing an overall factor of space-time volume, we 
define the effective potential Vcs {(j)c) as, 



d'^xVeBi(pc) 



Using now the definition of Dirac (5-function, 
and ([l5|) in (jl3|) we obtain, 



(16) 



(17) 



(18) 



Replacing (|18| ) in (14) we can write the effective action for constant field con- 
figurations as. 



r('^c) = £ ^C(27r)^5(^Ho)r(")(p. = 0) = f; ^cr(")(p, = o) /d^ 

ri=0 n=0 

and comparing it with ( p^ we obtain the final expression, 

Kff(</>e) = -^-Cr(")(K = o) 



a; 

(19) 
(20) 



Tl = 



which will be used for explicit calculations of the effective potential. 

Let us finally mention that there is an alternative way of expanding the 
effective action: it can also be expanded in powers of momentum, about the 
point where all external momenta vanish. In configuration space that expan- 
sion reads as: 



m = / d 



^^eff(0) + ^(a^0(a:))2^(0) + 



(21) 



1.2 The one-loop effective potential 

We are now ready to compute the effective potential. In particular the zero- 
loop contribution is simply the classical (tree-level) potential. The one-loop 
contribution is readily computed using the previous techniques and can be 
written in closed form for any field theory containing spinless particles, spin-^i 
fermions and gauge bosons. Here we will follow closely the calculation of Ref.ff. 
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Scalar fields 



We consider the simplest model of one self-interacting real scalar field, de- 
scribed by the lagrangian 



with a tree-level potential 



/: = ia^(/.9^0-yo(</') (22) 



V, = im^^^ -I- (23) 



The one-loop correction to the tree-level potential should be computed 
as the sum of all IPl diagrams with a single loop and zero external momenta. 
Diagrammatically they are displayed in Fig. |l|, where each vertex has 2 external 
legs. 



Figure 1: IPI diagrams contributing to the one-loop effective potential of (^). 

The n-th diagram has n propagators, n vertices and 2n external legs. The 
n propagators will contribute a factor of — m? + ie)~" |^. The external 
lines contribute a factor of 0^" and each vertex a factor of — iA/2, where the 
factor 1/2 comes from the fact that interchanging the 2 external lines of the 
vertex does not change the diagram. There is a global symmetry factor 
where ^ comes from the symmetry of the diagram under the discrete group of 
rotations Z„ and i from the symmetry of the diagram under reflection. Finally 
there is an integration over the loop momentum and an extra global factor of 
i from the definition of the generating functional. 

Using the previous rules the one-loop effective potential can be computed 

as, 

Kff(0c) = Vb((/)e)+l^l(0c), 



"We are using the Bjorken and Drell'sH notation and conventions. 



4 



with 



Vi{4>c) 



n=l 

i 

"2 



d'^p 1 

(27r)4 2^ 



A</>2/2 



(2^) 



■log 



p2 _ ^ ,jg 



After a Wick rotation 

/ = ^p% PE = {^^P^P), p" = ipy-p' = -pI, 



Eq. (gj) can be cast as, 



1 f d^PE 
(27r)4 



log- 



Finally, using the shifted mass 
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m {4>c) = m + -A0; 



d^Voi^c) 



(24) 



(25) 



(26) 



(27) 



and dropping the subindex E from the euclidean momenta, we can write the 
final expression of the one-loop contribution to the effective potential as. 



d'^p 

(2^ 



log[p^ +m\(f>^)] 



(28) 



where a field independent term has been neglected. 

The result of Eq. ( |2^ ) can be trivially generalized to the case of com- 
plex scalar fields described by the lagrangian. 



(29) 



The one-loop contribution to the effective potential in the theory described by 
the lagrangian (^^ is given by 



= l-Tr 



d'^p 
(27r)4 



log 



where 



(30) 



(31) 



and Tr = 2 V^, where the factor of 2 comes from the fact that each 
complex field contains two degrees of freedom. Similarly Tr 1=2 Ns- 



Fermion fields 

We consider now a theory with fermion fields described by the lagrangian, 

C = ilp,^-dr-^ai.Mf)l^'' (32) 
where the mass matrix {Mf)'^{<l)\) is a function of the scalar fields linear in 0^: 

The diagrams contributing to the one-loop effective potential are depicted 
in Fig. |. 






Figure 2: IPI diagrams contributing to the one- loop effective potential of (^)- 

Diagrams with an odd number of vertices are zero because of the 7- 



matrices property: tr{'-f^^ . . .-y 



0. The diagram with 2n vertices has 



2n fermionic propagators. The propagators yield a factor 

where Tr^ refers to spinor indices. The vertices contribute as 

rr[-i2»M/((?!)c)2"] 

where Tr runs over the different fermionic fields. There is also a combinatorial 
factor (from the cyclic and anticyclic symmetry of diagrams) and an overall 
— 1 coming from the fermions loop. One finally obtains the total factor 



1 rr(M^") 



2n 



■ Tr,l. 



The factor Tr^l just counts the number of degrees of freedom of the fermions. 
It is equal to 4 if Dirac fermions are used, and 2 if Weyl fermions (and a- 
matrices) are present. So we will write, 



Tr,l = 2A 



(33) 
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where A = 1 (A = 2) for Weyl (Dirac) fermfons. On the other hand we have 
grouped terms pairwise in the matrix product and used, 



-2 2 

p =p 



where p stands either for p • 7 or p • ct, depending on the kind of fermions we 
are using. 

Cohecting everything together we can write the one-loop contribution to 
the effective potential from fcrmion fields as, 



(27r)4 2^ 



Ml 

p2 



d^p 



(2n 



4 log 



Ml 

p2 



(34) 

As in the case of the scalar theory, after making a Wick rotation to the 
Euclidean momenta space, and neglecting an irrelevant field independent term, 
we can cast (BJ) as 



-2X-Tr 
2 



d'^p 



Iog[p2+M2(0,)] 



(35) 



Gauge bosons 

Consider now a theory described by the lagrangian, 



(36) 



In the Landau gauge, which does not require ghost-compensating terms, the 
free gauge-boson propagator is 



m, = -- 



p^ + ie 



with 



P^Pu 
p2 



satisfying the property p^A^'j, = and A" = A, n = 1, 2, 
The only vertex which contributes to one-loop is 



C 



1 



{M,,)l,AlA^P 



where 
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(37) 
(38) 

(39) 
(40) 



Figure 3: IPI diagrams contributing to the one-loop effective potential of (Bq). 



In this way the diagrams contributing to the one-loop effective potential are 
depicted in Fig. ||. 

A few comments about Eq. (^0|): (i) is the gauge coupling constant 
associated to the gauge field A^; if the gauge group is simple, e.g. SU{5), 
SO{lQ), Eq, . . ., then all gauge couplings are equal; otherwise there is a distinct 
gauge coupling per group factor, (ii) Ta are the generators of the Lie algebra 



of the gauge group in the representation of the (f>-&e\ds and the trace in (40) 
is over indices of that representation. 

Taking into account the combinatorial factors, the graph with n propaga- 
tors and n vertices yields a total factor 

2n p^" 

where 

Tr{A) = 3 (41) 

which is the number of degrees of freedom of a massive gauge boson. Collecting 
together all factors, and making the Wick rotation to the euclidean momenta 
space, we can cast the effective potential from gauge bosons as, 

V, = Tr(A)irr J ^ log [p' + (A/,,)2(0,)] (42) 



1.3 The higher-loop effective potential 

Calculating the effective potential by summing infinite series of Feynman graphs 
at zero external momentum is an extremely onerous task beyond the one-loop 
approximation. However, as has been shown in Ref. this task is trivial for 



the case of one-loop, and affordable for the case of higher-loop. Here we will 
just summarize the result of Ref.Q ^. 

We will start considering the theory described by a real scalar field, with a 
lagrangian Z given in (|2^p3|) , and an action as in ([|) . We will define another 
lagrangian £ by the following procedure: 



(43) 



where we have used in the last term the notation (||). In ([43|), (^c is an x- 
independent shifting field. The second term in (^) makes the vacuum energy 
equal to zero, and the third term is there to cancel the tadpole part of the 
shifted action. 

If we denote by 2?{0c; a; — v\ the propagator of the shifted theory. 



S(j){x)S(j){y) 



(44) 



and 



its Fourier transform, the effective potential is found to be given byE 



2 / (27r)4 H0c;p} 



I < exp 



i / d**x£/{0c;0(a;)} 



(45) 



The first term in ( [45|) is just the classical tree-level potential. The sec- 
ond term is the one-loop potential, where the determinant operates on any 
possible internal indices defining the propagator. The third term summarizes 
the following operation: Compute all IPI vacuum diagrams, with conventional 
Feynman rules, using the propagator of the shifted theory I?{(/>c;p} and the 
interaction provided by the interaction lagrangian £/{0c; 0(a;)}, and delete the 
overall factor of space-time volume [ d*x from the effective action (|l^). It can 
be shown that the last term in (^) starts at two- loop. Every term in (45) 
resums an infinite number of Feynman diagrams of the unshifted theory. 

In the simple example of the lagrangian ( p2[]2^ ) it can be easily seen that 
the shifted potential is given by 



1 



V{(bc; 0} - -m^(0,)0^ + -6,^' + -A' 



''The interested reader can find in0 all calculational details 



A 

4!' 



(46) 
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where the shifted mass is defined in (^7|). The shifted propagator is found to 
be 

iV-^{4>,-p}=p'' -m\4>,) (47) 
and the second term of ( p5[ ) easily computed to be 

Vi{<i^c) log[p' - ^'(<^c)] (48) 



which is easily seen to coincide with (^4|), up to field independent terms, so 
that after the Wick rotation we recover the result (p8|). 

The two-loop effective potential is harder than the one-loop term, but 
affordable. The result can be found in Ref.Q. Diagrammatically, the one- and 
two-loop effective potentials are given by the "figure-eight" plus the "sunset" 
diagrams, where it is understood that we are using the Feynman rules of the 
shifted theory, as stated above. 

Of course, the previous rules apply also to theories containing fermions 
and gauge bosons. The Feynman rules of the shifted theory applied to all IPI 
diagrams provide the total effective potential according to (^5|). In particular 
it is trivial to obtain the one-loop effective potential for fermions and gauge 
bosons, as given by Eqs. ( |35| ) and (42), respectively. Notice that the masses 
M|(0c), in (^), and M|^((/)c), in (42) are the masses in the corresponding 



shifted theories. Diagrammatically (35) and ( [4 21) can be represented in the 



shifted theory as vacuum diagrams with one fermion and one gauge boson 
loop, respectively. 



1.4 Renormalizations conditions 

The final expression of the effective potential we have deduced in the previous 
) and (^), is ultraviolet-divergent. To make sense out of 
it we have to follow the renormalization procedure of quantum field theories. 
First of all, to give a sense to the ultraviolet behaviour of the theory we have 
to make it finite: i.e. we have to regularize the theory. Second of all, all 
infinities have to be absorbed by appropriate counterterms, which were not 
explicitly written in our previous expressions. The way these infinities are 
absorbed by the counterterms depend on the definition of the renormalized 
parameters, i.e. on the choice of the renormalization conditions. Finally 
the theory, written as a function of the renormalized parameters, is finite. 

In this way, the first step towards renormalizing the theory is choosing the 
regularization scheme. We will first present the straightforward regularization 
using a cut-off of momenta. 



section, Eqs. (|28|), (35 
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Cut-ofT regularization 

We will illustrate this scheme with the simplest theory: a massless real scalar 
field, with a lagrangian 

£ = 1(1 + 6Z)id,^)' \5m^^- - (49) 

where 5Z , 5rn? and are the usual wave- function, mass and coupling constant 
renormalization counterterms. They have to be defined self-consistently order 
by order in the loop expansion. Here we will compute everything to one-loop 
order. 

The conventional definition of the renormalized mass of the scalar field is 
the negative inverse propagator at zero momentum. In view of ( pO| ) we can 
write it as: 

ml, = -r(2)(p = 0) = 5 (50) 

We can also define the renormalized coupling as the four-point function at zero 
external momentum, 



A« = -r(^)(p = o) = ^ 



(51) 

and the standard condition for the field renormalization is, 

Z(0) = 1 (52) 

Now we will compute the effective potential (|2^) cutting off the integral 
at = A^. First of all wc can integrate over angular variables. For that we 
can use, 

7r"/2 f 

dyf{p)^^J f{p)p''^'-'dp (53) 
where p — and we can cast (p8| ) as 

^i^^c) = ^l p\og[p + m'{<j,,)]dp. (54) 
327r^ Jo 

This indefinite integral can be solved with the help of 

1 1 fx^ 

xlog(a + x)dx — 2 (^^ ~ '^^) -'^'^S(fl + ^) ^ 2 ( "2 — 
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Neglecting now in ( ]5^ ) field independent terms, and terms which vanish 
in the hmit A — > oo, we finaUy obtain, 



1 



327r' 



rm2(0,)A2 



1 



647r2 



log 



A2 



(55) 



Using now (|55|) we can write the one-loop effective potential of the theory 
(H) as. 



,4 ^ 

' 64^2 



2567r2 



. A02 1 
^°^2A^-2 



(56) 



We will impose now a variant of the renormalization conditions (| 
and (b2). For the renormalized mass we can impose it to vanish, i.e., 







(51) 



(57) 



/>c=0 



For the renormalized gauge coupling A, we cannot use Eq. ( |51| ) at a value of 
the field equal to zero. There is nothing wrong with using a different renor- 
malization prescription and using a different subtraction point. We can use. 



A 



(58) 



where (x is some mass scale. Different choices of the scale lead to different 
definitions of the coupling constant, i.e. to different parametrizations of the 
same theory, but in principle any value of /i is as good as any other. 

Imposing now the conditions ( |57| ) and ( p8| ) to ( p6[ ) we can write the coun- 
terterms as, 

^ rA2 (59) 



and 



5\ = -- 



Sni^ = — 
11A2 



327r2 



327r2 

3A2 
32^ 



2A2 



(60) 



Using now ( p9| ) and (|60|) in (g6|) we can write the one-loop effective poten- 
tial in the previous renormalization scheme as. 



2567r2 



, , . 25 
log ^ - IT 



(61) 



A similar renormalization scheme can be defined also for theories with 
fermions and/or gauge bosons. However for gauge theories the regularization 
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provided by the cut-ofF explicitly break gauge invariance so that the dimen- 
sional regularization is better suited for them. In the next section we will 
review the calculation of the effective potential in the dimensional regulariza- 
tion and define the so-called MS scheme. 

Dimensional regularization 

This regularization scheme was introduced by t'Hooft and Veltmani. It con- 
sists in making an analytic continuation of Feynman integrals to the complex 
plane in the number of space-time dimensions n. The integrals have singu- 
larities which arise as poles in l/(n — 4) and have to be subtracted out. The 
particular prescription for subtraction is called a renormalization scheme. 
In working with the effective potential it is customary to use the so-called MS 
renormalization scheme 

We will compute now the one-loop effective potential (^Sj) using dimen- 
sional regularization, i.e. 



where ^ is a scale with mass dimension which needs to be introduced to balance 
the dimension of the integration measure. It is simpler to compute the tadpole 



where the meaning of V is the derivative with respect to m?{(t)c), using the 
basic formula of dimensional regularization, 



and integrating the resulting integral with respect to m?{(j)c). One can then 
write the regularized potential (|2|) as. 




(62) 




(63) 




(64) 




(65) 



We can expand ( |65|) in powers of 2 — rt/2 and use the expansion 

r(^) ^-~^E + o{z) 



(66) 
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where = 0.5772 ... is the Euler-Masccheroni constanti. We obtam for (65) 



647r2 



1 

2~~ 



7b + log 47r 



(67) 



Now the MS renormahzation scheme consists in subtracting the term pro 
portional to 

7b + fog 47r 



uv 



9 n 

^ 2 



(68) 



in the regularized potential (p7h. Therefore the divergent piece, 



647r2 



7_B + log Att 



has to be absorbed by the counterterms. Therefore the final expression for the 
one-loop potential, written in terms of the renormalized parameters, is 



1 



647r' 



:m'^{(j>c) < lo, 



9 



(69) 



For instance, in the theory described by lagrangian (^), the counterterms 
are given by, 



327r2 



9 n 

^ 2 



(5™^ = 
S\ = 

and the effective potential is, 

^^^=4! ^^+ 256^ 



(70) 



— JE + log 47r 



3 



(71) 



The scale /i along this section is related to the renormahzation group be- 
haviour of the renormalized couplings and masses. 

For a theory with fermion fields, one needs a trace operation in dimensional 
regularization, as Trl — f{n). For instance, for an even dimension one could 
choose, f{n) = 2"/^ foj- Dirac fermions, and f{n) — 2"/^-! for Weyl fermions. 
However the difference f{n) — /(4) is only relevant for divergent graphs and 
can therefore be absorbed by a renormalization-group transformation. It is 
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usually convenient to choose f{n) — /(4) = 2A for all values of 71 0. The 
effective potential (p5|) can be computed as in (E^), leading to, 



M|(0c) 



7_B + log 47r 



lOf 



(72) 



In the MS renormalization scheme, after subtracting the term proportional to 
(pq) we obtain. 



1 



^i(<^c) = -A^M;(0e)<;io: 



(73) 



Similarly, in a theory with gauge bosons as in (|3^), the effective potential 
(E3) is computed as, 



1 



o n 

^ 2 



- - + log 47r 



(74) 



where 



Tr{A) = n - 1 



(75) 



In the MS renormalization scheme, subtracting as usual the term propor- 
tional to (pq) one obtains the effective potential. 



= 3- 



1 



6471^ 



rM,\(0.) <^ lof 



(76) 



A variant of the MS renormalization scheme is the DR renormalization 
scheme til, where the dimensional regularization is applied only to the scalar 
part of the integrals, while all fermion and tensor indices are considered in four 
dimensions. In this case Tr{A) is taken equal to 3, as in (|4l|), and subtracting 
from (ItJ) the term proportional to (|68| ) one obtains. 



Vi( 



647r2 



(77) 
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1.5 One-loop effective potential for the Standard Model 

In this subsection we will apply the above ideas to compute the one loop 
effective potential for the Standard Model of electroweak interactions. The 
spin- zero fields of the Standard Model are described by the SU{2) doublet, 

Xi + iX2 

$ = I <l>c + h + ix3 I (78) 

where 0c is the real constant background, h the Higgs field, and Xa (a=li2,3) 
are the three Goldstone bosons. The tree level potential reads, in terms of the 
background field, as 

Vo{<Pc)^^^cbl + ^cbt (79) 
with positive A and , and the tree level minimum corresponding to 

The spin-zero field dependent masses are 

''^fii^'c) = 3A0C — 

m\{ct>,) = Ml~m^ (80) 

so that m\{v) = 2Aw^ = 2m^ and m^{v) = 0. The gauge bosons contributing 
to the one-loop effective potential are and Z, with tree level field dependent 
masses, 



i-W\ 



=) = ^'/'c (81) 

X _ .g' + .9" ,2 



Finally, the only fermion which can give a significant contribution to the one 
loop effective potential is the top quark, with a field-dependent mass 

mU<t>c) = ^cfi (82) 

where ht is the top quark Yukawa coupling. 

The one-loop effective potential Vi{4>c.) can be computed using Eqs. (p8|), 
(p5|) and (B2h. As we have said in the previous subsection, these integrals are 
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ultraviolet divergent. They have to be regularized and the divergent contribu- 
tions cancelled by the counterterms 

V,-'- =Sn + + ^-^c^t (83) 

where we have introduced a counterterm SQ for the vacuum energy or cosmo- 
logical constant (see next section). 

The final expression for the effective potential is finite and depends on the 
used regularization and, correspondingly, on the renormalization conditions. 
Next we will describe the two most commonly used renormalization conditions 
for the Standard Model. 



MS renormalization 



In this case we can use Eqs. (|67|), ( |72| ) and ( |74| ) for the contribution to Vi{(j)c) 
of the scalars, fermions and gauge bosons, respectively. In the MS renormaliza- 
tion scheme we subtract the terms proportional to Cuv , see Eq. ( |68| ) , which are 
cancelled by the counterterms in (|8^). One easily arrives to the finite effective 
potential provided by 



1 



647r2 



E 



i=W,ZJi,x,t 

where Ci are constants given by. 



Cw — Cz 



log 



mf (0c 



and Hi are the degrees of freedom 

nw =6, nz — 3, Uh = 1, n-^— 3, rit = —12 
The counterterms which cancel the infinities are provided by, 



(84) 



(85) 



(86) 



59. 
Srn? 



647r2 
3Ato2 



(nh + n^) Cuv 



167r2 



nh 



Ci 



uv 



167r2 



2g^ + {g^+ g'^f 
16 



(87) 



Cuv 
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where Cuv is defined in (|68|). We have explicitly written in the contribu- 
tion to the counterterms from the Higgs sector, and n^. The latter give rise 
entirely to the mass counterterms Sm^ and SQ. For Higgs masses lighter than 
W masses, the Higgs sector can be ignored in the one loop radiative corrections 
(as it is usually done) and the massive counterterms are not generated. 



Cut-ofT regularization 

A very useful scheme El is obtained by regularizing the theory with a cut-off 
and imposing that the minimum, at u = 246.22 GeV, and the Higgs mass does 
not change with respect to their tree level values, i.e., 



d^Vi+Vf- 



Now we can use (E55|) to write 
1 



327r^ 



E 



i=W,Z,tJi,x 



= 

= 



log 



A2 



(89) 



Imposing now the conditions (^) the infinities in ( |89| ) cancel against those 
in Vf'^', and the resulting 0c-dependent potential is finite, and given by, 



TTl j ((/)c) lo, 



mf(w) 



--\+2mi{v)mi{4>,) 



The counterterms (557, 5m^ and 5\ in (p3|) turn out to be given by 



5\ 
5rn? 



IGtt 



IGtt 



-y 

i 

-y 



mj{v) - hi 



log 



A2 



(90) 



(91) 



-mi{v)+b, log 



A2 



327r2 



E 



A'-m2(i;) + -fe, log 



where h 



w 



i=h,x 

bz = bt = and bh — b^ 



A2 



2 

-m . 
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We can see again in Q that ignoring the contribution to the one loop 
effective potential from the Higgs sector results in the non appearance of a 
cosmological constant. However, unlike the MS scheme, dm^ is also generated 
by the contribution of the gauge boson and top quark loops. Of course the one 
loop counterterms we are computing along this section are only useful for two 
loop calculations. 



1.6 Improved effective potential and renormalization group 



As we have seen in the previous section, the calculation of the effective action 
involves a mass /i which is not physical in the sense that all the theory should 
be independent of the chosen value of fi. In fact a change in /i should be ac- 
companied by a change in the renormalized parameters (couplings and masses) 
such that all the theory remains unchanged. This statement for the effective 



action can be expressed as an equationli 



^ dXi 







(92) 



for an appropriate choice of the coefficients Pi and 7, where denotes collec- 
tively all couplings and masses of the theory. In the last term of (|9^ ) we have 
made use of the notation (^) . 

We define the effective potential V as in Eq. (po|). 



00 

V = vifi, K, </.,) = i>(/i, A„ 0) - V -cr(") [p, = 0) 



(93) 



The role of the vacuum energy f2, 

f) - ^(m, A„ 0) 

has been recently stressed in Ref. El. Using now the renormalization group 



equation (RGE) satisfied by the effective action (^2[), we obtain the RGE sat- 
isfied by V as 



djjL ^ dXi d4>c 



V 



■ d 



f3^ 



d 



If we make a (/)c-independent shift to V such that, 

v = v + An{fi,x,) 
fi = h + An 



(94) 



(95) 
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with the condition, 



d ^ d - 



17 = 



then the potential V satisfies the well known RGE, 

V = {) 



d ^ d ,9 

Ofi OXi 0(j)c 



The formal solutions to Eqs. (^6|) and ( p7| ) can be written as, 

v = vi^l,x,,(f>c) = v{^l{t),x,it),<j>{t)) 



where 



n 


= Xi) = 


m 


= /iexp(t) 
= M{t) 


m 


- exp{- 


A(A(t)) 


dX^{t) 
dt 



/ ^{X,{t'))dt' 



with the boundary conditions, 



m = 1 

A,(0) - A 



In fact Eqs. ( |9q ) and (|9^) can be simply written as, 





dt 



— V = 
dt 

which state that fl and V are scale-independent. Of course the same 
to all derivatives of V, 
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which by virtue of (p9h satisfies 



^^"^ = atr-Q^v{ii{t),K{t):m) (103) 

The RGE satisfied by can be obtained from (p7|) and the property, 



Qn 



d<f)c 



It is given by, 



d 



V 



in) _ 



n'yV 



in) 



(104) 



which imphes that 1/^"^ is scale independent. 

In particular the scale independence of V^'^\ n — 0,1, . . ., means that we 
can fix the scale t at any value, even (/)c-dependent. Suppose we fix t by the 
arbitrary conditions, 

t = t{(p,) = \og{f{4,,)/^i} (105) 

m = m'pc))<pc 



Using ( |105| ) we can write the effective potential and its derivatives (102) 
as (/)c-functions, 

y(0,) = V [f{^,),x,{t{4>c)),m4>c))] (106) 

and 

\/(")(^,) ^ ^{t{cj^,)r^^v{^Ji{t),K{t), m 

Using Eq. (104) one can easily prove that0. 



(107) 



t=t(0c) 



(108) 



Fixing the scale is a matter of convention. Fixing the scale, as we have 
just described, as a function of 0^ {i.e. giving different scales for different 
values of the field) is usually done to optimize the validity of the perturbative 
expansion, i.e. minimizing the value of radiative corrections to the effective 
poteptial around the minimum of the field. A very interesting result obtained in 
Ref.113 is: The RGE improved effective potential exact up to (next-to-leading)^ 
log order ^ is obtained using the L-loop effective potential and the (L+l)-loop 
RGE P-functions. 

'^The convention is (next-to-leading)" =leading, i.e. L = 0. For L = 1 the potential is exact 
to next-to-leading log. 
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2 Field Theory at Finite Temperature 



The formalism used in conventional quantum field theory is suitable to describe 
observables {e.g. cross-sections) measured in empty space-time, as particle 
interactions in an accelerator. However, in the early stages of the universe, at 
high temperature, the environment had a non-negligible matter and radiation 
density, making the hypotheses of conventional field theories impracticable. 
For that reason, under those circumstances, the methods of conventional field 
theories are no longer in use, and should be replaced by others, closer to 
thermodynamics, where the background state is a thermal bath. This field 
has been called field theory at finite temperature and it is extremely useful to 
study all phenomena which happened in the payly universe: phase Ixaiisitions, 
inflationaoz.^ cosmology, ... Excellent articles cSEj, review articles E3'ta'E3 and 
textbooksEHl exist which discuss different aspects of these issues. In this section 
we will review the main methods which will be useful for the theory of phase 
transitions at finite temperature. 

2.1 Grand- canonical ensemble 

In this section we shall give some definitions borrowed from thermodynamics 
and statistical mechanics. The microcanonical ensemble is used to describe 
an isolated system with fixed energy E, particle number N and volume V. The 
canonical ensemble describes a system in contact with a heat reservoir at 
temperature T: the energy can be exchanged between them and T, N and V 
are fixed. Finally, in the grand canonical ensemble the system can exchange 
energy and particles with the reservoir: T, V and the chemical potentials are 
fixed. 

Consider now a dynamical system characterized by a hamiltonian ^ H and 
a set of conserved (mutually commuting) charges Qa- The equilibrium state of 
the system at rest in the large volume V is described by the grand-canonical 
density operator 



where $ = logTr exp {— J^a ^aQa ~ PH} is called the Massieu function (Leg- 
endre transform of the entropy), a a and [3 are Lagrange multipliers given by 
(3 = T~^, a A = —(3^1 A, T is the temperature and ^a are the chemical poten- 
tials. 

■^All operators will be considered in the Heisenberg picture. 




(109) 
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Using (109) one defines the grand canonical average of an arbitrary 
operator O, as 

(O) EE Tr{Op) (110) 

satisfying the property (1) = 1. 

In the following of this section we will always consider the case of zero 
chemical potential. It will be re-introduced when necessary. 

2.2 Generating functionals 

We will start considering the case of a real scalar field 4'{x), carrying no charges 
(fiA = 0), with hamiltonian H, i.e. 

(f>{x)^e''"(f>iO,x)e-''" (111) 

where the time x'^ — t is analytically continued to the complex plane. 

We define the thermal Green function as the grand canonical average of 
the ordered product of the n field operators 

^ (xi , . . . , Xji 

) = {Tc<l>{xi),...,<l>{xn)) (112) 

where the Tc ordering means that fields should be ordered along the path C 
in the complex i-plane. For instance the product of two fields is defined as, 

Tc<P{x)<j>{y) = 9c{x'> - y°)0(a;)0(y) + Ociy' - a:°)0(y)0(x) (113) 

If we parameterize G as t = z{t), where t is a real parameter, Tc ordering 
means standard ordering along t. Therefore the step and delta functions can 
be given as 0cit) = 6'(r), Scit) = (9z/9r)"^ S{t). 

The rules of the functional formalism can be applied as usual, with the 
prescription 5j{y) /Sj{x) = 6cix^ —y^)5^^^ (x—y), and the generating functional 
Z'^lj] for the full Green functions, 

^"bl^E^ / d'xi...d^x,,j{xi)...3{x,,)G^''Hxi,...,x^) (114) 

can also be written as, 

Zf[j] = (Tcexpl^tJ^d^xjixmx)^^ (115) 

which is normalized to Z^[0] — (1) — 1, as in (110), and where the integral 
along t is supposed to follow the path C in the complex plane. 
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Similarly, the generating functional for connected Green functions W^[j] 
is defined as Z^[j] = exp{iW^[j]}, and the generating functional for IPI Green 
functions r^[0], by the Legendre transformation, 

r^m = wf'ij] - / d'J-^pAj^^) (116) 

Jc '^J[^) 

where the current j(x) is eliminated in favor of the classical field ip^x) as <p{x) = 
5Wf^[j]/5j{x). It follows that ST'^[f]/d'^{x) = -j{x), and ^(x) = (0(x)) is the 
grand canonical average of the field (j){x). 
Symmetry violation is signaled by 



= (117) 

3=0 

for a value of the field different from zero. 

As in field theory at zero temperature, in a translationally invariant theory 
4>{x) = (j)c is SL constant. In this case, by removing the overall factor of space- 
time volume arising in each term of [0c] , we can define the effective potential 
at finite temperature as, 

T^m^~ Jd^xVfM (118) 
and symmetry breaking occurs when 



(119) 



for (Pc ^ 0. 



2.3 Green functions 
Scalar fields 

Not all the contours are allowed if we require Green functions to be analytic 



with respect to t. Using (113) we can write the two-point Green function as, 

G(^)(x -y) = edx' - y")G+{x - y) + dc{y" - xO)G_(x - y) (120) 
where 

G+ix~y)^{<Pix)<Piy)), G^x - y) ^ G+{y - x) (121) 
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Now, take the complete set of states \n) with eigenvalues i?„: H\n) 
En\n). One can readily compute (121) at the point x — y = as 



G+(x" 



y 



e-''^\{m\m\n)fe 



(122) 



m,n 



so that the convergence of the sum implies that —/3 < Im{x'^ — y") < which 
requires 9c{x^ — y°) = for Im{x^ — y^) > 0. From (121) it follows that the 
similar property for the convergence of G-{x^ —y^) is that < Im{x^ —y^) < /3, 
which requires 9c{y^ — x^) = for Im(iP — y^) < 0, and the final condition 
for the convergence of the complete Green function on the strip 

-P< Im{x° -y°) <P (123) 

is that we define the function 6c{t) such that 9c{t) — for Im{t) > 0. The 
latter condition implies that C must be such that a point moving along it has 
a monotonically decreasing or constant imaginary part. 

A very important periodicity relation affecting Green functions can be 
easily deduced from the very definition of G+{x) and G_(a;), Eq. (121). By 
using the definition of the grand canonical average and the cyclic permutation 
property of the trace of a product of operators, it can be easily deduced, 

G+{t~ip,x) ^ G-{t,x) (124) 

which is known as the Kubo-Martin-Schwinger relationE3. 

We can now compute the two-point Green function ( |120[ ) for a free scalar 
field, 

d^p 



4>{x) 



(2^)3/2(2Wp)l/2 



l{p)t 



a\p)t 



(125) 



where ujn 



y^p^ + TO^, which satisfies the equation 

7^2] G(^) {x-y) = -i5c{x - y) EE ~i5c{x'' - y'^)5^'^^ [x ~ y) (126) 



Using the time derivative of (|125|), and the equal time commutation rela- 



tion. 



4>{t,x),i{t,y) \ ^i6^^\x-y) 



(127) 



one easily obtains the commutation relation for creation and annihilation op- 
erators, 

[aip),aHk)]=S^'\p-k) (128) 
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and defining the Haniiltonian of the field as, 

H = l ^Lo,a\p)a(p) (129) 

one can obtain, using ( |128| ) the thermodynamical averages, 

{a\p)a{k)) - nB{uj.p)5^''\p- k) (130) 
{a{p)a^{k)) = [1 + nB{u:p)]5^''\p- k) 
where nB{Lo) is the Bose distribution function, 

(131) 



We will give here a simplified derivation of expression ( |13C| ). Consider 
the simpler example of a quantum mechanical state occupied by bosons of the 
same energy uj. There may be any number of bosons in that state and no 
interaction between the particles: we will denote that state by \n). The set 
{|n)} is complete. Creation and annihilation operators are denoted by and 
a, respectively. They act on the states |7i) as, a^\n) = \/rr+T|n + 1) and 
a\n) = \pn\n — 1), and satisfy the commutation relation, [a, a^] = 1. The 
hamiltonian and number operators are defined a.s H — loN and N — a^a, with 
eigenvalues ton and n, respectively. 

It is very easy to compute now (a^a) and (aa^) as in ( |130| ) using the 
completeness of {1?^)}. In particular. 



and 



n=0 n=0 



(l_e-/3-')2 



n=0 

from where (a^a) ~ nsioj), and (aa^) = 1 + nsii^), as we wanted to prove. 
Using now ( |130[ ) we can cast the two-point Green function as, 

G"^)(x - y) = / (0TP(p)e-'^("-^^ [^c(x° - y") + nsip")] (132) 

where the function p{p) is defined by p{p) — 27 r[0(p °) — 0(— p°)](5(p^ — m^). 
Now the particular value of the Green function ( |132| ) depends on the chosen 
contour C. We will show later on two particular contours giving rise to the 
so-called imaginary and real time formalisms. Before coming to them we will 
describe how the previous formulae apply to the case of fermion fields. 
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Fermion fields 

We will replace here (|l2^) and ( pH] ) by, 

s%\^ -y)^ {TcM^)^piy)) = edx' - y')s+^ - 0c(y° - x')s~^ (133) 

where a and /? are spinor indices, and 

S+,{x~y)^{Mx)i^fiiy)) (134) 

are the reduced Green function, which satisfy the Kubo-Martin-Schwinger re- 
lation, 

5^^(i-*/3,^)--^c:/3(i,^) (135) 
The calculation of the two-Green function for a free fermion field, satisfying 
the equation 

{ij-d~ m)^^ Sf^^ {x - y) = iScix - y)daf3 (136) 



follows lines similar to Eqs. (125) to ( |132| ). In particular, one can define a 
Green function S''^'^^ as 

5^5^ {x-y) = {ij-d + mUS^''^ {x ~ y) (137) 



where S'^'-^\x — y) satisfies the Klein-Gordon propagator equation (126). One 
can obtain for 5''^-' the expression, 

S^'^Kx -y) = j ^Mp^"^""'^ [^c(a;° - 2/°) - npip^)] (138) 
where 77.^(1^) is the Fermi distribution hmction 



Eq. ( |139| ) can be derived similarly to ( |131| ) as the mean number of fermions 
for a Fermi gas. This time the Pauli exclusion principle forbids more than one 
fermion occupying a single state, so that only the states |0) and |1) exist. They 
are acted on by creation and annihilation operators W and 6, respectively as: 
&'f|0) = |1), 5t|l) — 0, 5|0) = 0, 5|1) — |0), and satisfy anticommutation rules, 
{6, 6^} = 1. Defining the hamiltonian and number operators as H = luN and 
N — b^'b, we can compute now the statistical averages of (6^6) and (fefe^) using 
the completeness of {|n)}. 

n=0 n=0 
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and 



n=0 



from where (6^6) — np{u!), and = 1 — as we wanted to prove. 



2.4 Imaginary time formalism 

The calculation of the propagators in the previous sections depends on the 
chosen path C going from an initial arbitrary time t to t — i/3, provided by 
the Kubo-Martin-Schwinger periodicity properties (124) and (135) of Green 
functions. The simplest path is to take a straight line alonfithe imaginary axis 
t = —ir. It is called Matsubara contour, since Matsubarac3 was the first to set 
up a perturbation theory based upon this contour. In that case Sc{t) — ?(5(r). 

The two-point Green functions for scalar (132) and fermion ( |138| ) fields 
can be written as, 



G(r,f) 



(2^ 



p{p)e'P^e'^P° [e{T) + 7?n(/)] 



(140) 



where the symbol 77 stands as: rjs = ^ (rjp = —1) for bosons (for fermions). 
Analogously, n{p'^) stands either for nsip'''), as given by (131) for bosons, or 
npijp^), as given by (|139|) for fermions. It can be defined as a function of 77 as. 



(141) 



The Green function (140) can be decomposed as in (12C) 



G(T,f) = G+{T,x)e{T) + G^{T,x)e{-T) 



(142) 



Using now the Kubo-Martin-Schwinger relations, Eqs. (124) and (135), we 
can write G(r + (3) = r^Gir) for - /3 < t < 0, G(t - 13) ^ 7^G{t) for < 
T < [3, which means that the propagator for bosons (fermions) is periodic 
(antiperiodic) in the time variable r, with period [3. 

It follows that the Fourier transform of (140) 



G(w„,p)= / dr / d^xe*'^"^-"PG(T,x) 



(143) 



(where < a < /3) is independent of a and the discrete frequencies satisfy the 
relation rje^'^"^^ — 1, i.e. ujn ^ 2m: (3^^ for bosons, and Un = (2n + l)7r/3^^ for 
fermions. 
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Inserting now ( |140| ) into ( |143| ) we can obtain the propagator in momentum 
space G 

G{<^n,p) = 2 - (144) 

We can now define the euchdean propagator, A(— ir, a;), by 

G'(T,f) = iA(-ir,f) (145) 



where G{t,x) is the propagator defined in ( |140[ ). Therefore, using (144), we 
can write the inverse Fourier transformation, 

A{x) = i y / (146) 

P ^ J (27r)3 p2^^2_,_^2 V ' 

where the Mats ubara frequencies w„ are either for bosons or for fermions. 

From (146) one can deduce the Feynman rules for the different fields in 
the imaginary time formalism. We can summarize them in the following way: 

Boson propagator : -^r—^ — -; — \2nin(3~^ ,p] 

Fermion propagator : ; p'^ — \(2n + l)iTTP^^ ,p] 

J ■ p — m 

Loop integral : ^ / (147) 

n— — OO ^ / 

Vertex function : -iP{2T:)^Sj2^/^\J2P''> 



There is a standard trick to perform infinite summations as in (147). For 
the case of bosons we can have frequency sums as, 

oo 

- /(/=*^n) (148) 

n— — oo 

with w„ = 2n-Kj3~^ . Since the function -/3coth(-/3z) has poles at z = iujn and 
is analytic and bounded everywhere else, we can write (148) as, 

^ ^ dz/(z)^coth(i/3z) 



2mPJ^ '2 '2 
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where the contour 7 encircles anticlockwise all the previous poles of the imag- 
inary axis. We are assuming that /(z) does not have singularities along the 
imaginary axis (otherwise the previous expression is obviously not correct). 
The contour 7 can be deformed to a new contour consisting in two straight 
lines: the first one starting at —ioo + e and going to ioo + e, and the second 
one starting at ioo — e and ending at —ioo — e. Rearranging the exponentials 
in the hyperbolic cotangent one can write the previous expression as, 



1 
2^ 



dz^[f{z) + fi-z)] + ^^ 



dz[fiz) + fi-z)]- 



— -ioo+e 



,0z _ I 



and the contour of the second integral can be deformed to a contour C which 
encircles clockwise all singularities of the functions /(z) and f{—z) in the right 
half plane. Therefore we can write (148) as 



1 °° r 
- J2 f{p' = ^^n)^ j 



n— — oo 



ioo 1 r rl 

-^[/(z) + /(-z)]+ / -±nB{z)[f{z) + f{-z)] 

(149) 



where nsiz) is the Bose distribution function (131). 

Eq. (|149|) can be generalized for both bosons and fermions as, 



1 , 

- J2 /(/=^^n)= J ^ 



dz 



Ani 



.[/(z) + /(-z)]+77 



c 



dz 



n{z)[f{z) + f{-z)] 
(150) 



where the distribution functions n(z) are defined in (141). Eq. (15C) shows 
that the frequency sum naturally separates into a T independent piece, which 
should coincide with the similar quantity computed in the field theory at zero 
temperature, and a T dependent piece which vanishes in the limit T —^ 0, i.e. 

/3^oo. 



2. 5 Real time formalism 

The obvious disadvantage of the imaginary time formalism is to compute Green 
functions along imaginary time, so that going to the real time has to be done 
through a process of analytic continuation. However, a direct evaluation of 
Green function in the real time is possible by a judicious choice of the contour 
C in (112). The family of such real time contours is depicted in Fig. |^ where 



the contour C is C = Ci |J C2 |J C3 |J C4 where Ci goes from the initial time 
ti to the final time t/, C3 from tj to tj — ia, with < cr < /3, C2 from 
tf — ia to ti — ia, and C4 from U — ia to U — ij3. Different choices of a lead 
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Figure 4: Contour used in the real time formalism. 

to an equivalence class of quantum field theories at finite temperatur/^^^' 
instance the choice a = leads to the Keldysh perturbation expansioni 
the choice a = P/2 is the preferred one to compute Green functions. 




Computing the Green function for scalar (132) and fermion (138) fields 
taking path C is a matter of calculation, as we did for the imaginary time for- 
malism in ( |l40| )-( p^4[ ). One can, nrove that the contribution from the contours 
C3 and C4 can be neglected OJEa. Therefore, for the propagator between 
and y'^ there are four possibilities depending on whether they are on Ci or C2 ■ 
Correspondingly, there are four propagators which are labeled by (11), (12), 
(21) and (22). 

Making the choice a = (3/2, the propagators for scalar fields (132) can be 
written, in momentum space, as 



G(p) 



G(ii)(p) G(i2)(p) 

G(21)(p) G(22)(p) 



A{p) 






A*{p) 



Mb{P,p) 



(151) 

where A(p) is the boson propagator at zero temperature, and the matrix 
Mb{P,p) is given by. 



cosh 9 {p) smh9{p) 
sinh6'(p) cosh 9 (p) 



(152) 



where 



sinh6'(p) 



(1 - e-^"-) 
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cosh9{p) = (l — e 



(153) 



Using now (151), ( |152| ), ( |153D , one can easily write the expression for the 
four bosonic propagators, as 



^(21) ^ g(12) 



(154) 



Similarly, the propagators for fermion fields can be written as 

S{jp)al3 = 

= Mf{P,p), q 
where the matrix Mp{P,p) is, 
Mf{P,p) 



(7 • p + m)ai3/^ip) 

{'yp + m)af}A*{p) 



with 



cos6'(p) sin6'(p) 
sin6'(p) cos6'(p) 



cos dip) = [0(pO) - 0(-p°)] (1 + e-'3-p) -1/2 



(155) 

MfW.p) 

(156) 



(157) 



In the same way, using now ( |155| ), (156), and (157) one can easily write 
the expression for the four fermionic propagators, as 



2 2 
m 



)) 



5(")(p) = {^■p + m){A{p)-27TnF{u;p)5{p^ 

5(22) ^ 5(11)* (158) 

^(12) ^ ~2Tr{j ■ p + m)[0{p'') - 9i-p°)]e^'^''/^nF{ujp)6ip^ -m^) 

^(21) ^ _gil2) 



As one can see from (154) and (15S), the main feature of the real time 
formalism is that the propagators come in two terms: one which is the same 
as in the zero temperature field theory, and a second one where all the tem- 
perature dependence is contained. This is welcome. However the propagators 
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(12), (21) and (22) arc unphysical since one of their time arguments has an 
imaginary component. They are required for the cons isten cy of the theory. 
The only physical propagator is the (11) component in (|l54| ) and (158). 

Now the Feynman rules in the real time formalism are very similar to those 
in the zero temperature field theory. In fact all diagrams have the same topol- 
ogy as in the zero temperature field theory and the same symmetry factors. 
However, associated to every field there are two possible vertices, 1 and 2, and 
four possible propagators, (11), (12), (21) and (22) connecting them. All of 
them have to be considered for the consistency of the theory. In the Feynman 
rules, type 2 vertices are hermitian conjugate with respect to type 1 vertices. 
The golden rule is that: physical legs must always be attached to type 1 vertices. 
Apart from the previous prescription, one must sum over all the configurations 
of type 1 and type 2 vertices, and use the propagator G*^"^' or S'^"''^ to connect 
vertex a with vertex h. 

There is now a general agreement in the sense that the imaginary tipifi 
formalism and the real time formalism should give the same physical answercj. 
Using one or the other is sometimes a matter of taste, though in some cases 
the choice is dictated by calculational simplicity depending on the physical 
problem one is dealing with. 



3 The effective potential at finite temperature 

In this section we will construct the (one-loop) effective potential at finite 
temperature, using all the tools provided in the previous sections. As we will 
see, in particular, the effective potential at finite temperature contains the 
effective potential at zero temperature computed in section 1. The usefulness 
of this construction is addressed to the theory of phase transitions at finite 
temperature. The latter being essential for the understanding of phenomena 
as: inflation, baryon asymmetry generation, quark-gluon plasma transition 
in QCD,... We will compare different methods leading to the same result, 
including the use of both the imaginary and the real time formalisms. This 
exercise can be useful mainly to face more complicated problems than those 
which will be developed in this course. 



3. 1 Scalar fields 



We will consider here the simplest model of one self-interacting scalar fields 
described by the lagrangian (^2|) and (|23|). We have to co mput e the diagrams 
contained in Fig. ^ using the Feynman rules described in (147), for the imag- 
inary time formalism, or in (154) for the real time formalism. We will write 
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the result as, 

^off 

where Vo{(j)c) is the tree level potential 



lfff(0c) = Vo(0c) + lf ('^c) (159) 



Imaginary time formalism 

We will compute the diagrams in Fig. |^. Using the Feynman rules in Eq. ( |l47| ), 
Eq. (^) translates into, 

71— — OO V / 

where w„ are the bosonic Matsubara frequencies and 

uj^ = p ^ + {(l)c) (161) 

being defined in (p7|). 



The sum over n in ( |16C| ) diverges, but the infinite part does not depend 
on (j)c. The finite part, mtich contains the (/)c dependence, can be computed 
by the following method ta. Define, 

oo 

v{iu)^ log(^^+c.2) (162) 

then, 



n— — oo 



E (163) 

Tl— — oo 



Using the identity. 



oo ^ ^ 

fiy) = J2^ir^ = -l^ + l^^<^oth7ry (164) 



1 TT e-^'^y 



- + TT 



2y 2 1 - e-^'^y 



with y ~ (310/211 we obtain. 



dv 



1 e-P" 

2 ^ 1 - e-/^'^ 
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(165) 



and 



v{uj) = 2/3 



- + -log(l-e-^-) 



to ~ independent terms 



Substituting finally (|166|) into (16C) one gets, 



Vfi<f>c) 



cP'p 
(2^ 



| + ilog(l-e-) 



(166) 



(167) 



One can easily prove that the first integral in (167) is the one-loop effective 
potential at zero temperature. For that we have to prove the identity, 

/ — log(-a;'^ + - «e) = - + constant (168) 

2 J-oo 27r 2 



dx 



1 



(169) 



Integral (|169|) can be performed closing the integration interval (— cx3,oo) in 
the complex x plane along a contour going anticlockwise and picking the pole 
of the integrand at a; = —^JllP' — ie with a residue 1/2cl;. Using the residues 
theorem Eq. (16£) can be easily checked. Now we can use identity (168) to 
write the temperature independent part of (167) as 



(27r)3 



(2^)4 



log(-p2 + _ -^^ 



and, after making the Wick rotation p'^ = ipE in (170) we obtain. 



d^p _ 1 
(2^'^ ~ 2 



(2^ 



log[p2+m2 (,/),)] 



(170) 



(171) 



which is the same result we obtained in the zero temperature field theory, see 
Now the temperature dependent part in (167) can be easily written as, 



1 



d3p 



log (1 - e 



1 



/3 J (27r)3 ^ ' 27r2/34 

where the thermal bosonic function Jb is defined as. 



JB[m'(0c)/3' 



2 o2i 



JB[m'l3 



dx X log 



1 - e 



(172) 



(173) 
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The integral ( |173[ ) and therefore the thermal bosonic effective potential 
admits a high-temperature expansion which will be very useful for practical 
applications. It is given by 



A 2 2 / 9 \ 3/2 4 



' 45 ^ 12 r2 6 V T2 j 32 abT"^ 



(174) 



,C(2^ + 2) 
(^+1)! 



1 \ / \ 



2 y \^t:^T^ ) 



where a?, = 167r2exp(3/2 — 27^) (log at = 5.4076) and C, is the Riemann C- 
function. 

There is a very simple way of computing the effective potential: it consists 
in computing its derivative in the shifted theory and then integrating^. In 
fact the derivative of the effective potential 

dV,^ 



is described d iagra mmatically by the tadpole diagram. In fact using the Feyn- 
man rules in (147) one can easily write for the tadpole the expression, 



dyf__\(j)cl ^ f d^p 1 

J^^J (2^)3c.2+c.2 



n— — 00 



(175) 



or, using the expression (|27| ) for m'^{(l)c), 



dVf 



dm?{<i3c) 2f3 



^ OG 



71 — — 00 



d^p 



(27r)3 w2 + ^2 



(176) 



Now we can perform the infinite sum in (lit) using the result in Eq. (149) 
with a function / defined as, 



1 



(177) 



and obtain for the tadpole ( |176| ) the result 
dVf f d^p fl dz 1 



dm'^{(j)c) J (27r)3 [2 J 2T:i lu^ - J(j 2t:i e^^■' - L uj^ ~ z" 



dz 



1 



27ri e^^ - 1 ^2 



(178) 
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The first term in (178) gives the /3-independent part of the tadpole contribution 
as, 

1 r°° dz 1 

2 ./__ 2^cj2_z2 



(179) 



We can now close the integration contour of ( |l 79| ) anticlockwise and pick the 
pole of (177) a.t z — —uj with a residue l/2a;. The result of (|179|) is 



4a; 



(180) 



The second term in (178) gives the /3-dependent part of the tadpole contribu- 
tion. Here the integration contour encircles the pole a,t z — uj with a residue 



1 



1 



2w ef^"^ 



Adding (|180|) and (181) we obtain for the tadpole the final expression. 



dm'^{(j>c) 



d^p 
(2^ 



1 

2^ 



(181) 



(182) 



Now, integration of ( 182 ) with respect to m?{(t)c) leads to the expression (167) 
for the thermal effective potential and, therefore, to the final expression given 
by (p7l|) and (|l7|). 



Real time formalism 



As we will see in this section, the final result for the effective potential (167) 
can be also obtained using the real time formalism. Let us compute the tadpole 
diagram. Since physical legs must be attached to type 1 vertices, the vertex 
in the tadpole must be considered of type 1, and the propagator circulating 
around the loop has to be considered as a (11) propagator. Application of the 
Feynman rules (154) to the tadpole diagram leads to the expression H 



dV,^ 



A0C 
2 



(27r)4 



27mB(w)(5(p^ - w?{4>c)) 



(183) 



or, using as before the expression (|27|) for m {(f>c) 



dV,^ 



d-m?{4}c) 



d'^p 



(184) 



■^We are replacing in (154) the value of ujp by the corresponding value u given by (161) in 
the shifted theory. 
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Now the /3-independent part of (184), after integration on 



tributes to the effective potential as 



(2^ 



log(-p2 + m'^{(j)c) - ie) 



con- 



(185) 



Finally using Eq. ( |168| ) to perform the p° integral, we can cast Eq. (185) as 



(27r)3 2 



(186) 



which coincides with the first term in (167). 

Integration over p^ in the /3-dependent part of (184) can be easily per- 
formed with the help of the identity 



1 



)^ — [6{p°+iu,) + S{p'~Lu,)] 



leading to, 



d^p 1 



(187) 



(188) 



(27r)3 2io 

which, upon integration over m'^{(j)c) leads to the second term of Eq. ( |167| ). 

We have checked that trivially the real time and imaginary time formalisms 
lead to the same expression of the thermal effective potential, in the one loop 
approximation. 



3. 2 Fermion fields 

We will consider here a theory with fermion fields described by the lagrangian 
(^). As in the scalar case, we have to compute the diagrams contained in 
Fig. ^, using the Feynman rules either for the imaginary or for the real time 
formalism, and decompose the thermal effective potential as in (159). 



Imaginary time formalism 



The calculation of the diagrams in Fig. g, using the Feynman rules (147), 
yields. 



2A 

'2/3 



E 



n— — oo 



d^p 
(2^ 



where cun are the fermionic Matsubara frequencies and 



Ml 



(189) 



(190) 
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The sum over n is done with the help of the same trick employed in ( \16z 
(|166|). Let f{y) be given by (164), then, 



E 



y 



m=2,4,... 



E 



m— 1,3,... 



and using (164) we get, 



E 



y 



^ w2 + 4n2 2 

n=l " 
/(.) - \f (f ) 



TT TT 1 

4 2 e-^v + 1 



m— 1,3 

The function v{uj) in this case can be written as 



v(uj)^2 log 



/?2 



(191) 



(192) 



(193) 



and its derivative. 



dv _ 4/3 ^ y 
du) TT y^ + V? 



where y = Poj/t:. Then using (192) we get 



dv 
duj 



2(3 



1 



1 + e^'^ 



and, after integration with respect to w. 



v{uj) = 2/3 



Replacing finally (|196| ) into (|189| ) one gets. 



-2A 



(27r)3 



w 1 
2+^ 



log (1 + e 



(194) 



(195) 



a; — independent terms (196) 



(197) 



1681)- (p7l|), to lead to the one- 



The first integral in (197) can be proven, as in 

loop effective potential at zero temperature (|35). The second integral, which 
contains all the temperature dependent part, can be written as, 



2Ai 



log (1 + e-^'^) 



-2A- 



1 



■27r2/34 



(198) 
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where the thermal fcrmionic function Jp is defined as, 

1 + 



dx X log 



(199) 



As in the scalar field, the integral ( |199| ) and therefore the thermal fermionic 
effective potential admits a high-temperature expansion which will be very 
useful for practical applications. It is given by 



t\ I I 360 24 r2 32 T4 

r7/2 



(200) 

2 \ ^+2 
2^2 



where ay = 7r^exp(3/2 — 2^e) (logo/ — 2.6351) and C is the Riemann C,- 
function. 

As we did in the case of the scalar field, there is a very simple way of 
obtaining the effective potential, computing the tadpole in the shifted theory, 
and integrating over 0c. Using for the fermion propagator (147) 

.7 • p + Mf 
* p2 _ M] 

and the trace formula, 

Tri-y-p + Mf) = 2AM/ 
we can write for the tadpole the expression. 



n— — oo 



d^p 



(27r)3 Lul + uj^ 



(201) 



or, using the expression AIf{(j3c) — T^c, where T is the Yukawa coupling. 



dV- 



dM]{ 



-2A- 



1 °° 



2/3 



n— — oo 



d'^p 



(27r)3 w2 + c^2 



(202) 



Now the infinite sum in ( |202| ) can be done with the help of (15C), with 
f{z) given by (pT7|), as 



dVf 



dIvPMc) 



-2A 



(27r)3 [2 



1 



27ri a;2 — 2:2 



1 



1 



27ri e^'^ + 1 ^2 _ ^2 

(203) 
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The first term of (p03[) reproduces the zero temperature result (|35|), after 



Mf integration, by closing the integration contour of (179) anticlockwise and 



picking the pole at z = —lu with a residue l/2u;. The second term in (203) 
gives the /3-dependent part of the tadpole contribution. Here the integration 
contour C encircles the pole at z = uj with a residue 



(-2A) 



1 



1 



2UJ 6/3" + 1 

Adding all of them together, we obtain for the tadpole the final expression 



(204) 



= -A 



(2^ 



1 1 



1 



2uj uj eP'^ + 1 



(205) 



and, upon integration with respect to My we obtain the result previously 
presented in Eq. (197). 



Real time formalism 



As for the case of scalar fields, the thermal effective potential for fermions (197) 
can also be very easily obtained using the real time formalism. We compute 
again the tadpole diagram, where the vertex between the two fermions and 
the scalar is of type 1 and the fermion propagator circulating along the loop 
is a (11) propagator. Apphcation of the Feynman rules (158) leads to the 
expression 



-TTr 



(2^ 



(7-p + My) 



- Mj + ie 

or, using as before the expression for M?, 



27mF(w)(5(r -M/) 



(206) 



dVl' 



Trl 



dMfi 



d^p 



- 2TmFiio)6{p'^ - Mf) 



(207) 



Now the /3-independent part of (207), after integration on My, contributes 
to the effective potential, 

d^p U! 



Trl 



(27r)3 2 



(208) 



which coincides with the first term in (197). 
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Integration over in the /3-dependent part of (207) can be easily per- 



formed with the help of the identity (187) leading to 



which, upon integration over M| leads to the second term of Eq. (197) 



3.3 Gauge bosons 

The thermal effective potential for gauge bosons in a theory described by the 
lagrangian ( ^6| ) is computed in the same way as for previous fields. The sim- 
plest thing is to compute the tadpole diagram using the shifted mass for the 
gauge boson. In the Landau gauge, the gauge boson propagator reads as, 

nM^(p)("/3) ^ a^^g("^)(p) (210) 



where A is the projector defined in ( pSD with a trace equal to 3 (see Eq. (|4lD). 
Therefore the final expression for the thermal effective potential is computed 
as, 

yf(</),)=Tr(A)|i J -0^log[p' + M'^M] + ^JbIM^MP']] 

(211) 



where the thermal bosonic function Jb in ( |173| ). The first term of (211) agrees 
with the zero temperature effective potential computed in ( ^2| ) , and the second 
one just counts that of a scalar field theory a number of times equal to the 
number of degrees of freedom (3) of the gauge boson. 

3.4 The Standard Model case 

The Standard Model of electroweak interactions was previously defined through 
Eqs. (^)-(^), and the corresponding one loop effective potential at zero 
temperature computed through Eqs. (83)-(^l|) using various rcnormalization 



schemes and the contribution of gauge and Higgs bosons and the top quark 
fermion to radiative corrections. Here we will compute the corresponding one 
loop effective potential at finite temperature. We will use the rcnormalization 
scheme of Eq. (|8^), so that the renormalized effective potential at zero tem- 
perature is given by Eq. ( |90| ) , and consider only the contribution of W and Z 
bosons, and the top quark to radiative corrections. This is expected to be a 
good enough approximation for Higgs masses lighter than the W mass. 
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Using Eqs. ( |198| ) and (211) one can easily see that the finite-temperature 
part of the one- loop effective potential can be written as, 



2^ 



n,JB[mU4>c)/T^]+ntJF[niUc^c)/T^ 



=w,z 



(212) 



where the function Jb and Jp where defined in Eqs. (173) and ( |199| ), respec- 
tively. 

Using now the high temperature expansions ( |174D and ( |200| ), and the one 
loop effective potential at zero temperature, Eq. (pQ), one can write the total 
potential as. 



V{<Pc, T) = D{T^ ~ T^)4>1 ~ ET4)l 



A(T) 



(213) 



where the coefficients are given by 



D = 



2mi 



2mi 



8^2 



m c 



B 



^ O 

3 



il - 8Bv^ 



{2m 



4D 

4 I 
W + 



Am 



A(T) = A- 



(214) 

(215) 

(216) 
(217) 



2m^Y log 



-I- TO 7 log 



Am^ log 



9 



(218) 

where lo g Ab — lo g afc — 3/2 and logAp — loga^ — 3/2, and as, ap are 
given in ( 174 ) and(200). All t he masses which appear in the definition of the 
coefficients, Eqs. ( |214| ) to ( 218 ), are the physical masses, i.e. the masses at the 
zero temperature minimum. The peculiar form of the potential, as given by 
Eq. ( ^13[ ), will be useful to study the associated phase transition, as we will 
see in subsequent sections. 



4 Cosmological phase transitions 

All cosmological applications of field theories are based on the theory of phase 
transitions at finite temperature, that we will briefly describe throughout this 
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section. The main point here is that at finite temperature, the equiHbrium 
value of the scalar field (fi, ((/)(T)), does not correspond to the minimum of 
the effective potential V^'^'^{(j)), but to the minimum of the finite temperature 
effective potential Vfg{(f>), as given by (15£). Thus, even if the minimum of 
Vgg-^°((/()) occurs at (0) = cr 7^ 0, very often, for sufficiently large temperatures, 
the minimum of 1/^^(0) occurs at {(f>{T)) = 0: this phenomenon is known 
as symmetry restoration at high temperature, and gives rise to the nhase 
transition from (j){T) = to = u. It was discovered by Kirzhnits EZl in 
the context of the electroweak theory (symmetry breaking between weak and 
electromagnetic interactions occurs when the universe cools down to a critical 
temperature XL.r^. 1.0^. GeV) and subsequently confirmed and developed by 
other authorsB'Byy. 

The cosmological scenario can be drawn as follows: In the theory of the 
hot big bang, the universe is initially at very high temperature and, depending 
on the function Vfg{(j)), it can be in the symmetric phase {4>{T)) = 0, i.e. 
= can be the stable absolute minimum. At some critical temperature 
Tc the minimum at = becomes metastable and the phase transition may 
proceed. The phase transition may be first or second order. First-order phase 
transitions have supercooled (out of equilibrium) symmetric states when the 
temperature decreases and are of use for baryogenesis purposes. Second-order 
phase transitions are used in the so-called new inflationary modelscSI. We will 
illustrate these kinds of phase transitions with very simple examples. 

4.1 First and second order phase transitions 

We will illustrate the difference between first and second order phase transi- 
tions by considering first the simple example of a potential ^ described by the 
function, 

F(0, T) = D{T^ - t2)02 + (219) 
where D and are constant terms and A is a slowly varying function of T ^ 



A quick glance at (174) and (200) shows that the potential ( pi9[ ) can be part 
of the one- loop finite temperature effective potential in field theories. 

At zero temperature, the potential has a negative mass-squared term, 
which indicates that the state ^ = is unstable, and the energetically favored 

state corresponds to the minimum at 0(0) = ^^T^,, where the symmetry 
<f> <-> —(j) of the original theory is spontaneously broken. 



■''The fj) independent terms in (pia), i.e. ^^(0, T), are not explicitly considered. 
3 The T dependence of A will often be neglected in this section. 
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The curvature of the finite temperature potential ( |219| ) is now T-dependent, 
n? (0, T) ^ 3X^^ + 2D{T^ - T^) (220) 
and its stationary points, i.e. solutions to dV{(j),T) / d(j) — 0, given by, 

m = 

and (221) 



m = 



/2D(r2 -r2) 



\{T) 



Therefore the critical temperature is given by To. AtT > To, m^{0,T) > 
and the origin (/> = is a minimum. At the same time only the solution = 



in (221) does exist. At T = To, m'^{0,To) = and both solutions in {U2V 



collapse at (f) — 0. The potential ( |219| ) becomes, 

y(0,r„) = ^(/.4 (222) 
At T < To, m2(0, r) < and the origin becomes a maximum. Simultaneously, 



the solution 4>{T) ^ does appear in (221). This phase transition is called of 



second order, because there is no barrier between the symmetric and broken 
phases. Actually, when the broken phase is formed, the origin (symmetric 
phase) becomes a maximum. The phase transition may be achieved by a 
thermal fluctuation for a field located at the origin. 

However, in many interesting theories there is a barrier between the sym- 
metric and broken phases. This is characteristic of first order phase transi- 
tions. A typical example is provided by the potential ^, 

V{(l), T) = D{T^ - T^)(l)^ - ET4>^ + ^j^^'' (223) 
where, as before, D, Tq and E are T independent coefhcients, and A is a slowly 



varying T-dependent function. Notice that the difference between (223) and 



(219) is the cubic term with coefficient E. This term ca n be provided by the 
co ntrib ution to the effective potential of bosonic fields ( |l 74 |) The behaviour 
of ( ^23|) for the different temperatures is reviewed in Refs.li30. At T > Ti the 
only minimum is at (/) = 0. At T = Ti 

T^ ^^(T^^^^o (224) 
' ~ 8XiT,)D-9E^ ^ > 



''See, e.g. the one-loop effective potential for the Standard Model, Eq. (|90[) 
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a local minimum at 4>{T) ^ appears as an inflection point. The value of the 
field at T = Ti is, 

WTO) = ^ (225) 

A barrier between the latter and the minimum at = starts to develop at 
lower temperatures. Then the point ( ^25[ ) splits into a maximum 

SET 1 



(T) = 2^ - 2X(r) V^E^T^^WtJW^^) (226) 



and a local minimum 



SET 1 



0m (^) = ^A(T) + 2A(T) \/9^^'^' - 8A(T)Z?(T2 - T^) (227) 
At a given temperature T = Tc 

2 _ HTc)DT^ , . 

^ XiT,)D-E^ ^^^^^ 



the origin and the minimum (227) become degenerate. From (226) and (227) 
we find that 

ET 

MT.) = ^ (229) 

and 

0m(7^c) = ^ (230) 

For T <Tc the minimum at = becomes metastable and the minimum 
at (pmiT) ^ becomes the global one. At T = To the barrier disappears, the 
origin becomes a maximum 

^m{To)=0 (231) 
and the second minimum becomes equal to 

0m(ro) = |^ (232) 

The phase transition starts at T = Tc by tunneling. However, if the barrier is 
high enough the tunneling effect is very small and the phase transition does 
effectively start at a temperature Tj. > Tt > To- In some models To can be 
equal to zero. The details of the phase transition depend therefore on the 
process of tunneling from the false to the global minimum. These details will 
be studied in the rest of this section. 
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4-2 Thermal tunneling 

The transition from the false to the true vacuum proceeds via thermal tun- 
neling at finite temperature. It can be understood in terms of formation of 
bubbles of the broken phase in the sea of the symmetric phase. Once this has 
happened, the bubble spreads throughout the universe converting false vacuum 
into true one. ___ 

The tunneling jmteESSa is computed by using the rules of field theory at 
finite temperature Ej. In the previous section we defined the critical tempera- 
ture Tc as the temperature at which the two minima of the potential V{4>, T) 
have the same depth. However, tunneling with formation of bubbles of the 
field (p corresponding to the second minimum starts somewhat later, and goes 
sufficiently fast to fill the universe with bubbles of the new phase only at some 
lower temperature Tt when the corresponding eucMean action Se — S3/T 
suppressing the tunnehng becomes 0(130 — 140)EjO'E£|, as we will see in the 
next section. 

We will use as prototype the potential of Eq. ( ^23| ) which can trigger, 
as we showed in this section, a first order phase transition. In this case the 
false minimum is </) = 0, and the value of the potential at the origin is zero, 
l/(fl_T) = 0. The tunneling probability per unit time per unit volume is given 
hyM 

- - A{T)e-^''/^ (233) 

V 

In ( |233| ) the prefactor A(T) is roughly of 0(T'*) while ^3 is the three-dimensional 
euclidean action defined as 



5, = 



i(w)' + y(^,T) 



(234) 



At very high temperature the bounce solution has 0(3) symmetryEl and 
the euclidean action is then simplified to. 
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47r / r dr 
'q 



2 I dr 



V{P{r),T) 



(235) 



where — x and the euclidean equation of motion yields. 



with the boundary conditions 



dr'^ r ar 



lim 0(r) 

r — 'OO 



= 



(236) 
(237) 
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dr 



(238) 



r=0 



From here on we will follow the discussion in Ref. tHi. Let us take (f> — 
outside a bubble. Then ( |235| ), which is also the surplus free energy of a true 
vacuum bubble, can be written as 



R 



r^dr 



1 

2 V 



+ v{4>{t),t) 



(239) 



where R is the bubble radius. There are two contributions to (239): a surface 
term Fs, coming from the derivative term in ( ^39| ) , and a volume term _FV, 
coming from the second term in (p39h. They scale like, 



5, - 2ttR^ 



5j>_ 
SR 



6R 



AnR^iV) 



(240) 



and {V) is the average 



where SR is the thickness of the bubble wall, Stp — 
of the potential inside the bubble. 

For temperatures just below Tc, the height of the barrier V{(t)M, T) is large 
compared to the depth of the potential at the minimum, — y(0„i,T). In that 
case, the solution of minimal action corresponds to minimizing the contribution 
to Fy coming from the region (j) = (j)M- This amounts to a very small bubble 
wall SR/R <C 1 and so a very quick change of the field from (p = outside the 
bubble to (f) = 0^ inside the bubble. Therefore, the first formed bubbles after 
Tc are thin vifall bubbles. 

Subsequently, when the temperature drops towards To the height of the 
barrier V{(l)]\i, T) becomes small as compared with the depth of the potential at 
the minimum ~'V{(j)rm T). In that case the contribution to Fy from the region 
= is negligible, and the minimal action corresponds to minimizing the 
surface term Fg. This amounts to a configuration where 5R is as large as 
possible, i.e. 6R/R = 0{1): thick wall bubbles. So whether the phase 
transition proceeds through thin or thick wall bubbles depends on how large 
the bubble nucleation rate (|233| ) is, or how small 53 is, before thick bubbles 
are energetically favoured. i_. 

For the potential ( ^23| ) an analytic formula has been obtained in Ref. l3 
without assuming the thin wall approximation. It is given by. 



S3 
T 



13.72 



D 1- 



rp2 



3/2 



\{T)D 



1 - 



2^2 



(241) 
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1 



2.4 



0.26 



(242) 



l~x (l-a;)2_ 

The case of two fields is extremely more complicated. In particular the two- 
Higgs situation in the supersymmetric standard theory has been recently solved 
in Ref. E3. The connection between zero temperature and finite temperature 
tunneling is manifest. In particular at temperatures much less than the inverse 
radius the [T — 0) 0(4) solution has the least action. This can happen for 
theories with a supercooled symmetric phase: for instance in the presence of 
a barrier that does not disappear when the temperature drops to zero. At 
temperatures much larger than the inverse radius, the 0(3) solution has the 
least action. 



^.3 Bubble nucleation 

In the previous subsection we have established the free energy and the critical 
radius of a bubble large enough to grow after formation. The subsequent 
progress of the phase transition depends on th e ratio of the rate of production 



of bubbles of true vacuum, as given by ( 233 ), over the expansion rate of the 
universe. For example if the former remains always smaller than the latter, 
then the state will be trapped in the supercooled false vacuum. Otherwise the 
phase transition will start at some temperature Tt by bubble nucleation. 
The probability of bubble formation per unit time per unit volume is given by 



(HI) where B{T) = Sz(T)lT, A{T) = loT'^, where the parameter uj wiU be 
taken of 0(1). 

Since the progress of the phase transition should depend on the expansion 
rate of the universe, we have to describe the universe at temperatures close to 
the electroweak phase transition. A homogeneous and isotropic (flat) universe 
is described by a Robertson- Walker metric which, in comoving coordinates, is 
given by ds^ = (M?' — a(i)^ (dr^ -f r^dfi^) ^ where a{t) is the scale factor of the 
universe. The universe expansion is governed by the equation 



Mill 



P (243) 



where Afp^ is the Planck mass, and p is the energy density. For temperatures 
T ^ 10^ GeV the universe is radiation dominated, and its energy density is 
given by, 



. 30 (244) 

where g{T) = gsiT) + lgF{T), and gsiT) [gpiT)) is the effective number of 
bosonic (fermionic) degrees of freedom at the temperature T. For the stan- 
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dard model we have g^^^ = 106.75 which can be considered as temperature 

independent. 

The equation of motion ( ^43|) can be solved, and assuming an adiabatic 
expansion of the universe, a{Ti)Ti — a{T2)T2, one obtains the following rela- 
tionship, 

t = (245) 

where C = -;— a/ — ~ 3 x 10^^. Using (|245|) the horizon volume is given by 
47r V T^g 

VH{t)^K'^ (246) 

The onset of nucleation happens at a temperature Tt such that the prob- 
ability for_a. single bubble to be nucleated within one horizon volume 
is ~ 1, i.e. 123 

[ f (^)'exp{-^3(T)m = 0[l) . (247) 

which implies numerically, 

B(r.,.,37 + ,„gi^ + 41ogl52|£^ (248) 

where we have normalized Tt ~ 100 GeV and E^/{XD) ~ 10~^ which are 
typical values which will be obtained in the standard model of clectroweak 
interactions. 



5 Baryogenesis at phase transitions 

There are two essential problems to be understood related with the baryon 
number of the universe: 

i) There is no evidence of antimatter in the universe. In fact, there is no 
antimatter in the solar system, and only p in cosmic rays. However antiprotons 
can be produced as secondaries in collisions (pp —^^p + p) at a rate similar 
to the observed one. Numerically, — ^ 3 x lO^'', and — — ^ 10^^. We can 

conclude that 3> n-g, so tt-as = ns — n-g ^ ns- 

ii) The second problem is to understand the origin of 

7?= — - (0.3- 1.0) X 10-^ (249) 
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today. This parameter is essential for primordial nucleosynthesis l2L 77 may not 
have changed since nucleosynthesis. At these energy scales (~ 1 MeV) baryon 
number is conserved if there are no processes which would have produced 
entropy to change the photon number. We can easily estimate from rj the 
baryon to entropy ratio by using 



^3.91 = 7.04 



(250) 



and the range (249). 

In the standard cosmological model there is no explanation for the small- 
ness of the ratio (249) if we start from tiab = 0. An initial asymmetry has 
to be imposed by hand as an initial condition (which violates any naturalness 
principle) or has to be dynamically generated at phase transitions, which is 
the way we will explore all along this section. 



5.1 Conditions for baryogenesis 

As we have seen in the previous subsection the universe was initially baryon 
symmetric (jib — n-^) although the matter- antimatter asymmetry appears to 
be large today (uab ~ ns 3> n-^). In the standard cosmological model there is 
no explanation for the value of -q consistent with nucleosynthesis, Eq. ( |249| ), and 
it has to be imposed byJiand as an initial condition. However, it was suggested 
by Sakharov long ago l2I that a tiny riAs might have been produced in the 
early universe leading, after pp annihilations, to (249). The three ingredients 
necessary for baryogenesis are: 



B-nonconserving interactions 

This condition is obvious since we want to start with a baryon symmetric 
universe (A_B — 0) and evolve it to a universe where AS 7^ 0. B-nonconserving 
interactions might mediate proton decay; in that case the phmomenological 
constraints are provided by the proton lifetime measurements Ej Tp > \Q^^yr. 



C and CP violation 

The action of C (charge conjugation) and CP (combined action of charge 
conjugation and parity) interchanges particles with antiparticles, changing 
therefore the sign of B. For instance if we describe spin-i fermions by two- 
component fields of definite chirality (left-handed fields ijjL and right-handed 
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fields V'/j) the action of C and CP over them is given by 



p 


- 




IpR > 




c 


IpL - 








CP 






V'i?^ — * 





If the universe is initiaUy matter- antimatter symmetric, and without a 
preferred direction of time as in the standard cosmological model, it is repre- 
sented by a C and CP invariant state, \(j)o), with B — 0. If C and CP were 
conserved, i.e. [C, H] = [CP, H] = Q, H being the hamiltonian, then the state 
of the universe at a later time t, \(j>{t)) = e*^*|(/)o) would be C and CP invariant 
and, therefore, baryon number conserving, AB = 0. The only way to generate 
a net AB 7^ is to have C and CP violating interactions. 



Departure from thermal equilibrium 

If all particles in the universe remained in thermal equilibrium, then no direc- 
tion for time would be defined and CPT invariance would prevent the appear- 
ance of any baryon excess, rendering CP violating interactions irrelevant e3. 

A particle species is in thermal equilibrium if all its reaction rates, F, are 
much faster than the expansion rate of the universe, H. On the other hand 
a departure from thermal equilibrium is expected whenever a rate crucial for 
maintaining it is less than the expansion rate (T < H) . Deviation from thermal 
equilibrium cannot occur in a homogeneous isotropic universe containing only 
massless species: massive species are needed in general for such deviations to 
occur. 



5.2 Bary agenesis at the electroweak phase transitions 

It has been recently realized @@ that the three Sakharov's conditions for 
baryogenesis can be fulfilled at the electroweak phase transition: 

• Baryonic charge non-conservation was discovered by 't HooftEl. In fact 
baryon and lepton number are conserved anomalous global symmetries 
in the Standard Model. They are violated by non-perturbative effects. 

• CP violation can be generated in the Standard Model from phases in the 
fermion mass matrix, Cabibbo, Kobayashi, Maskawa (CKM) phases l3. 
This effect is much too small to explain the observed baryon to entropy 
ratio. However, in extensions of the Standard Model as the minimal 
supersymmetric standard model (MSSM), a sizeable CP violation can 
happen through an extended Higgs sector. 
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• The out of equilibrium condition can be achieved, if the phase transition 
is strong enough first order, in the bubble walls. In that case the B- 
violating interactions are out of equilibrium in the bubble jsialls and a 
net B-number can be generated during the phase transitioned. 



Baryon and lepton number violation in the electroweak theory 

Violation of baryon and lepton number in the electroweak theory is a very 
striking phenomenon. Classically, baryonic and leptonic currents are conserved 
in the electroweak theory. However, that conservation is spoiled by quantum 
corrections through the chiral anomaly associated with triangle fermionic loop 
in external gauge fields. The calculation gives, 

9^^'B = d,f,=N,[£-^WW-£,YY) (252) 

where Nj is the number of fermion generations, W^u and Y^i, are the gauge 
field strength tensors for SU{2) and J7(l)y, respectively, and the tilde means 
the dual tensor. 



A very important feature of (252) is that the difference B — L is strictly 
conserved, and so only the sum B + L is anomalous and can be violated. 
Another feature is that fluctuations of the gauge field strengths can lead to 
fluctuations of the corresponding value oi B + L. The product of gauge field 



strengths on the right hand side of Eq. (252) can be written as four-divergences 



WW = d^k'^,YY^ a^fc^, where 

fc^ = e^'^'^^Y.^Yp (253) 

and Wf^, Y^ are the gauge fields of SU{2) and U{1)y, respectively. In general 
total derivatives are unobservable because they can be integrated by parts 
and drop from the integrals. This is true for the terms in the four-vectors 



(253) proportional to the field strengths W^^ and y^,y. This means that for 
the abelian subgroup U{1)y the current non conservation induced by quantum 
effects becomes non observable. However this is not mandatory for gauge fields, 
for which the integral can be nonzero. Hence only for non-abelian groups can 
the current non conservation induced by quantum effects become observable. 
In particular one can write AB — AL = NfANcs, where Ncs is the so- 
called Chern-Simons number characterizing the topology of the gauge field 
configuration, 

2 /> 

Ncs^^J d'^e^'" {wtjW^ - ^-eabcW^W^wi) (254) 
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Note that though Ncs is not gauge invariant, its variation ANcs is- 

We want to compute now AB between an initial and a final configura- 
tion of gauge fields. We are considering vacuum field strength tensors 
which vanish. The corresponding potentials are not necessarily zero but can 
be represented by purely gauge fields, 

W^ = --U{x)d^U'\x) (255) 
9 

There are two classes of gauge transformations keeping W^^ = 0: 

• Continuous transformations of the potentials yielding ANcs = 0. 

• If one tries to generate ANcs 7^ by a continuous variation of the 
potentials, then one has to enter a region where W^^ 7^ 0. This means 
that vacuum states with different topological charges are separated by 
potential barriers. 

The probability of harrier penetration can be calculated using the quasi- 
classical approximation E3. In euclidean space time, the trajectory in field 
space configuration which connects two vacua differing by a unit of topological 
charge is called instanton. The euclidean action evaluated at this trajectory 
gives the probability for barrier penetration as F ~ exp (j--^^ ^ 10^^®^, 

where aw = 5^/47r. This number is so small that the calculation of the 
pre-exponential is unnecessary and the probability for barrier penetration is 
essentially zero. 



Baryon violation at finite temperature: sphalerons 

However, in a system with non zero temperature a particle may classically go 
over the barrier with a probability determined by the Boltzmann exponent, as 
we have seen. 

What we have is a potential which depends on the gauge field configuration 
Wfj_. This potential has an infinite number of degenerate minima, labeled as 
rin ■ These minima arc characterized by different values of the Chern-Simons 
number. The minimum Qq corresponds to the configuration = and we 
can take conventionally the value of the potential at this point to be zero. Other 



minima have gauge fields given by (255). In the temporal gauge Wq = 0, the 



gauge transformation U must be time independent (since we are considering 
gauge configurations with W/^^, — 0), i.e. U — U{x), and so functions U define 
maps, 

U -.S^ — > SU{2) 
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All the minima with W^^i, = have equally zero potential energy, but those 
defined by a map U{x) with nonzero Chern-Simons number 

n[U] = j d^xe'^''Tr{Ud,U-^UdjU-^UdkU-^) (256) 

correspond to degenerate minima in the configuration space with non-zero 
baryon and lepton number. 

Degenerate minima are separated by a potential barrier. The field config- 
uration at the top of the barrier is called sphaLeron, which is a static unstable 
solution to the classic equationS|-C|f motion c3. The sphaleron solution has 
been explicitly computed in Ref. eZI for the case of zero Weinberg angle, (i.e. 
neglecting terms 0{g')), and for an arbitrary value of sin^ 9w in Ref. ea. 

An ansatz for the sphaleron solution for the case of|-zero Weinberg angle 
was given (for the zero temperature potential) in Ref. 113, for the Standard 
Model with a single Higgs doublet, as, 

Wta'^dx' = -jfiOdU U-^ (257) 

for the gauge field, and 

$ = ^h{£,)U \ ) (258) 
for the Higgs field, where the gauge transformation U is taken to be, 

r \ —X + ly z J 

and we have introduced the dimensionless radial distance ^ = gvr. 

Using the ansatz (257), (258) and (25£) the field equations reduce to, 

= 2/(l-/)(l-2/)- jft'd-/) (260) 

|(«^|) - 2Mi-/)^ + A,^(.= -i). 

with the boundary conditions, /(O) = h{Q) = and f{oo) = h{oo) = 1. The 
energy functional becomes then, 
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The solution to Eqs. (260) has to be found numerically. The solutions 
depend on the gau ge a nd quartic couplings, g and A. Once replaced into the 
energy functional ( |26l| ) they give the sphaleron energy which is the height of 
the barrier between different degenerate minima. It is customary to write the 
solution as, 

i?.ph = ^S(A/g^) (262) 
aw 

where B is the constant which requires numerical evaluation. For the standard 
model with a single Higgs doublet this parameter ranges from B{{)) — 1.5 to 
B{oo) — 2.7. A fit vahd for values of the Higgs mass 25 GeV < mt < 250 GeV 
can be written as, 

B(a;) = 1.58 + 0.32x- 0.05x^ (263) 

where x — rrih/mw- 

The previous calculation of the sphaleron energy was performed at, ^rp 
temperature. The sphaleron at finite temperature was computed in Ref. mrH. 
Its energy follows the approximate scaling law, Esph{T) — £'sph(</'(2^))/(</>(0)) 



which, using (262), can be written as. 



iJ3ph(T) = ^H^BiX/g') (264) 
aw 

where mw{T) = ig(<^(T)) 
Baryon violation rate at T > 

The calculation of the baryon violation rate at T > Tc, i.e. in the sym- 
metric phase, is very different from that in the broken phase, that will be 
reviewed in the next section. In the symmetric phase, at </> = 0, the pertur- 
bation theory is spoiled by infrared divergences, and so we cannot rely upon 
perturbative calculations to compute the baryon violation rate in this phase. 
In fact, the infrared divergences are cut off by the non-perturbative genera- 
tion of a magnetic mass, niM ^ ctwT, i.e. a magnetic screening length 
{<^wT)^^. The rate of baryon violation per unit time and unit volume 
r does not contain any exponential Boltzmaiin factor |. The pre-exponential 
can be computed from dimensional groundsEj as 

r = k{awTf (265) 

where the coefficient k has been evaluated numerically in Ref. with the 
result 00.1 <k< 1.0. 



'It would disappear from ( p33| ) in thip4jmit T ^ oo. 

^In fact recent lattice computations H seem to provide an extra factor in (265) as k'a-w 
which is roughly speaking 0(1). 
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Baryon violation rate at T < Tc 

After the phase transition, the calculation of baryon violation rate can be done 



using the semiclassical approximations given by Eq. (233). The rate per unit 
time and unit volume for fluctuations between neighboring minima contains 
a Boltzmann suppression factor exp {—Esph{T)/T), where Esph{T) is given by 
(264), and a pre- factor containing the detepninant of all zero and non-zero 
modes. The prefactor was computed in Ref.Eil as 



2.8 X lO^T^ f ^ 
V 47r 



4 C 

K 



(266) 



where we have defined C,{T) — Esph{T)/T, the coefficient B is the function 
of X/g^ defined in ( ^63| ) and k is the functional determiaant associated with 
fluctuations about the sphaleron. It has been estimatedl2Zl to be in the range, 
10-4 < K< 10-1. 

The equation describing the dilutipa S of the baryon asymmetry in the 
anomalous electroweak processes readsEJ 



dS_ 
'dt 



-VB{t)S 



(267) 



where Vsit) is the rate of the baryon number non-conserving processes. As- 
suming T is constant during the phase transition the integration of (267) yields 
13., r 



S 



and X 



-N. 



T3 



t. Using now (266) and (245) we can write the 



exponent X as, X ^ lO^^K^^e-'', where we have taken the values of the pa- 
rameters, B — 1.87, aw = 0.0336, Nf = 3, Tc 10^ GeV . Imposing now the 
condition S ^ 10"^, or A" < 10, leads to the condition on C{Tc), 



C{T,) > 71ogC(rc) + 91ogl0 + logK 



(268) 



Now;, taking k at its upper bound, k — 10 ^, we obtain from ( ^68[ ) the 
bound E3 



Esph{Tc) ^ 



45, 



and using the lower bound, k = 10 we obtain, 

Esph{Tc) 



Tc 



37, 



(269) 



(270) 



Eq. (269) is the usual bound used to test different theories while Eq. (270) 
gives an idea on how much can one move away from the bound (269), i.e. the 
uncertainty on the bound (269). 
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The bounds ( |269D a nd ( |270| ) can be translated mto bounds on (f>{Tc)/Tc. 
Using the relation ( p64[ ) we can write 



4>iT,) _ g g,ph(T,) 1 g,ph(rc) 



36 



(271) 



where we have used the previous values of the parameters. The bound (269) 
translates into 

m) > 



while the bound ( ^70D translates into, 

m) > 



1.3 



1.0 



(272) 



(273) 



These bounds, Eqs. (272) and (269), require that the phase transition is 
strong enough first order. In fact for a second order phase transition, (j){Tc) ~ 
and any previously generated baryon asymmetry would be washed out during 
the phase transition. For the case of the Standard Model the previous bounds 
translate into a bound on the Higgs mass, as we will see. 



6 On the validity of the perturbative expansion 

The approach of Ref. to the finite temperature effective potential relied 
on the observation that symmetry restoration implies that ordinary 
perturbation theory must break down at high temperature. In fact, 
otherwise perturbation theory should hold and, since the tree level potential is 
temperature independent, radiative corrections (which are temperature depen- 
dent) should be unable to restore the symmetry. We will see that the failure 
of perturbative expansion is intimately linked to the appearance of infrared 
divergences for the zero Matsubara modes of bosonic degrees of freedom. This 
just means that the usual perturbative expansion in powers of the coupling 
constant fails at temperatures beyond the critical temperature. It has to be 
replaced by an improved perturbative expansion where an infinite number of 
diagrams are resummed at each order in the new expansion. We will review 
the actual situation in this section. 



6.1 The breakdown of perturbative expansion 

We will examine the simplest model of one self-interacting real scalar field, 
described by a lagrangian with a squared-mass term, rri^ and a quartic cou- 
pling A. We will use now power counting arguments to investigate the high 
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temperat 
potential '^'51 



-lehaviour of higher loop diagrams contributing to the effective 
Elca. After rescaling all loop momenta and energies by T, a loop 
amplitude with superficial divergence D takes the form, T^f{m/T). If there 
are no infrared divergences when m/T — > 0, then the loop goes like T^. For 
instance the diagram contributing to the self-energy of Fig. ^ is quadratically 




Figure 5: One-loop contribution to the self-energy for the scalar theory. 



divergent {D — 2), and so behaves like XT^. For D < 0, there are infrared 
divergences associated to the zero modes of bosonic propagators in the imagi- 
nary time formalism [n — in (147)] and the only T dependence comes from 
the T in front of the loop integral in (147). Then every logarithmically diver- 
gent or convergent loop contributes a factor of T. For instance the diagram 
contributing to the self-energy in Fig. contains two logarithmically divergent 




Figure 6: Two-loop contribution to the self-energy for the scalar theory. 

loops and so behaves like, X^T^ = A(Ar^). It is clear that to a fixed order in 
the loop expansion the largest graphs are those with the maximum number of 
quadratically divergent loops. These diagrams are obtained from the diagram 
in Fig. ^ by adding n quadratically divergent loops on top of it, as shown in 
Fig. I 

They behave as, 



(274) 



where M is the mass scale of the theory, and has been introduced to rescale 
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Figure 7: Daisy (n+l)-loop contribution to the self-energy for the scalar theory. 



the powers of the temperature ^ As was clear from Eq. (274), adding a 



quadratically divergent bubble to a propagator which is part of a logarithmi- 
cally divergent or finite loop amounts to multiplying the diagram by 

(275) 

This means that for the one-loop approximation to be valid it is required that 

A— < 1 

AP 

along with the usual requirement for the ordinary perturbation expansion 

A < 1 



However at the critical temperature we have that M/%/A [see Eqs. (pT9|)- 
( p20| )]. Therefore at the critical temperature the one-loop approx- 
imation is not vaUd and higher loop diagrams where multiple quadrati- 
calljj—divergent bubbles are inserted cannot be neglected. Daisy resumma- 
tionS consists precisely to resum all powers of a and provides a theory where 
mP'ii'c) fnctf = "^^(0c) +n, where 11 is the self-energy corresponding to the 
one-loop resummed diagrams to leading order in powers of the temperature 
T (e.g. ^ T^). This method was pursued systematically by ParwaniEj and 
applied to the Standard Model by Arnold and EspinosaEJ. 

What about the diagrams which are not considered in the improved ex- 
pansion? The two-loop diagram of Fig. ^ is suppressed with respect to the 
diagram of Fig. ^ by A. On the other hand the multiple loop diagram obtained 

*In fact, the mass M has a different meaning for the improved and the unimproved theories, 
as we shall see. For the unimproved theory, M is the mass in the shifted lagrangian, = 
m-^{4>), while in the improved theory, M is given by the thermal mass. 
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Figure 8: Non-daisy (n+2)-loop contribution to the self-energy for the scalar theory. 

from that of Fig. ^ by adding n quadratically divergent loops on top of it, see 
Fig. behaves as 



rp2n+2 rp2n+l / rp 

M2» M2»-i I M 



(276) 



and it is suppressed with respect to the multiple loop diagram of Eq. ( P74D 
by \T/M []. Therefore the validity of the improved expansion is guaranteed 
provided that, 



A< 1 



(277) 



This simplified discussion can be done of course in more complicated field 
theories, as e.g. the Standard Model, with similar arguments EJ. We will now 
apply these results to two different models. 



6.2 The scalar theory 

In the case just considered of the scalar theory, the effective potential can be 
written in the one-loop daisy resummed approximation as 

y(</>, T) = im,\02 _ J_ (^^2^ + 1 ^^2^ ' + ^ . . . (278) 
where m^g = + ^T^. 

'Non-daisy contributions to the self-e nerg y are suppressed, with respect to daisy contribu- 
tions, by O{0), where (3 is defined in (E77|). The correspondme contributions to the vacuum 
diagrams (i.e. effective potential) are suppressed by C'(/3^)Ea. 
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The potential (278) apparently yields a first order phase transition El. 
However, the symmetry breaking minimum occurs when the three terms in 
the potential are similar. Then, at T ^ Tc, at the minimum 



VXT, 



'cff 



and then, 



A- 



rrics 



which shows that the result of perturbation theory is fake 



E. 



(279) 
(280) 



6.3 The Standard Model 

In cases where there are gauge and/or Yukawa interactions, as the Standard 
Model case, the situation is completely differ ent. The effective potential in the 
one-loop improved approximation, Eq. (213) can be written as 



1 



A 



F(0, T) = -mis^' - bg'<j>''T + -^^ + 



(281) 



where 



bg'4,'T s [4M,=,(*) + 2MU4>)] -p^ 



9 



T 

12^ 



and rn^a 



ag^T^ , ag^ denoting the contribution of gauge and Yukawa 



couplings to the one loop self-energy. We have only considered the contribution 
of transverse gauge bosons to the phase transition strength, and neglected that 
of the (screened) longitudinal gauge bosons. 

Now again, symmetry breaking occurs when all terms are similar, 



9'^ rji 2 9^ rr2 

-T, ™.,~-r 



and the expansion parameter /3, 



T 



Mw{(l)) (t) g 



T A 
9-- — 



(282) 



(283) 



Therefore the validity of perturbation theory implies tiiat niH ~ raw ■ This 
behaviour is supported by non-perturbative calculationsE3. 
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7 Summary of physical results for the phase transition and con- 
strains on the Higgs mass 

7.1 Standard Model results 

The effective potential for the Standard Model was analyzed in Eqs. ( |213|) , 
(281) in the one-loop approximation, including leading order plasma effects. 
In this approximation, the longitudinal components of the gauge bosons are 
screened by plasma effects while the transverse components remain unscreened. 
In this way a good approximation to the effective potential including these 
plasma effects is provided by Eq. (223), where the coefficient E is given by 



22m|. 



9.5 X 10-^ 



(284) 



Now we can use Eq. (230) and m\ — 2\v^ to write, 

Tr. ^ ml 



(285) 



In this way the bound (272) translates into the bound on the Higgs mass. 



rrih 



' — - 42 GeV. 



(286) 



The bound (28£) is excluded by LEP measurementsEil, and so the Standard 
Model is unable to keep any previously generated baryon asynimetry. Including 
two- loop effects the bound is slightly increased to ~ 45 GeVlEl. Is it possible, 
in extensions of the Standard Model, to overcome this difficulty? We will see 
in the next section one typical example where the Standard Model is extended: 
the well motivated supersymmetric extension of the Standard Model. 



7.2 MSSM results 

Among the extensions of the Standard Model, the physically most motivated 
and phenomenologically most acceptable one is the Minimal Supersymmetric 
Standard Model. This model allows for extra CP-violating phases besides the 
Kobayashi-Maskawa one, which could help in generating the observed baryon 
asymmetry til. It is then interesting to study whether in the MSSM the na- 
ture of the phase transition can be significantly modified with respect to the 
Standard Model. 

In this section we extend the considerations of the previous sectioBu-to 
the full parameter space, characterizing the Higgs sector of the MSSMEa'E3e3. 
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The main tool for our study is the one-loop, daisy-improved finite-temperature 

effective potential of the MSSM, Vcs{4'jT). We are actually interested in the 
dependence of the potential on (pi = Re and (j)2 = Re H2 only, where 
and H2 are the neutral components of the Higgs doublets Hi and H2, thus 
(f> will stand for {(j)i,(f)2)- Working in the 't Hooft-Landau gauge and in the 
Di^-scheme, we can write 



Kff(</',T) 
where 



Vo (</>) + Vi 0) + AVi (</., T) + AKiaisy (0, T) + V2 {<P, T) , (287) 



Vo{<l>) 
Vi{<i>,Q) 

AVi{cP,T) 

Aydaisy(</',T) 



= mf 0f + m2?!>2 + 2rnl(j}i(j)2 + ■ 



,2 i2 




127r 



ni[Ml{cP,T)-ri4{4>)\ 



,(288) 

(289) 

(290) 
(291) 



The masses Mi{(l),T) are obtained from the m?(0) by adding the leading T- 
dependent self-energy contributions, which are proportional to T^. We recall 
that, in the gauge boson sector, only the longitudinal components {Wl, -^z,, 7i,) 
receive such contributions. 

The relevant degrees of freedom for our calculation are: rit = —12, nj^ = 
rif^ = 6, nw = 6, nz = 3, nw^ = 2, nz^^ = n-^^ = 1. The ficld-dcpcndcnt top 
mass is m1{4') = hlcf)^- The entries of the field-dependent stop mass matrix 
are 



.m?(</.) + D|(0), 



ht{At(f>2 - n4>i) 



(292) 
(293) 
(294) 



where mq, mu and At are soft supersymmetry-breaking mass parameters, /x 
is a superpotential Higgs mass term, and 



1 2 , 
- - - sm^ 

2 3 



■ sm U\Y 



I 2 



.9' + g 



5 Wl 



(0? - "/-a), 



(295) 
(296) 
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are the D-term contributions. The field-dependent stop masses are then 



77^ ((/.)+ ml (0) 



m~ {(l))-m^ {(/)) 



+ [m|(0)]^ (297) 



The corresponding effective T-dependent masses, Al~ {4>,T), are given by 

tl,2 

expressions identical to (297), apart from the replacement 



tL,R 



m4 



(298) 



The leading parts of the T-dependent self-energies of tL,R , 

n~(r) = lg2j,2^ 1^2 j,2^_l_ ^, 2^,2^1 /^2j.2^ (299) 
n~(r) = Ig2j,2^^^,2y2^ 1^2j,2^ (30Q^ 

where gs is the strong gauge coupling constant. Only loops of gauge bosons, 
Higgs bosons and third generation squarks have been included, implicitly as- 
suming that all remaining supersymmetric particles are heavy and decouple. 

We shall work in the limit in which the left handed stop is heavy, niQ ^ 500 
GeV. In this limit, the supersymmetric corrections to the precision electroweak 
parameter Ap become small and hence, this allows a good fit to the electroweak 
precision data coming from LEP and SLD. Lower values of niQ make the phase 
transition stronger and we are hence taking a conservative assumption from the 
point of view of defining the region consistent with electroweak baryogenesis. 
The left handed stop decouples at finite temperature, but, at zero temperature, 
it sets the scale of the Higgs masses as a function of tan/3. For right-handed 
stop masses below, or of order of, the top quark mass, and for large values 
of the CP-odd Higgs mass, tua ^ Mz, the one-loop improved Higgs effective 
potential admits a high temperature expansion, 

Vo+Vi = — -T 

(301) 

where Nc = 3 is the number of colours and i?sM is the cubic term coefficient 
in the Standard Model case. 



/ \a/2- 

12^ 



A(T) 



65 



Within our approximation, the lightest stop mass is approximately given 



by 



~ + 0.15M^ cos 213 + mill 



(302) 



where 
in Ref. 

breaking parameter = 



= At — n/ tern (3 is the stop mixing parameter. As was observed 
the phase transition strength is maximized for values of the soft 
Ilj:^{T), for which the coefficient of the cubic 
term in the effective potential, 



E~Ei 



SM 



3/2 



4y27r 



(303) 



can be one order of magnitude larger than -Esm Ea. In principle, the above 
would allow a sufhciently strong first order phase transition for Higgs masses 
as large as 100 GeV. However, it was also noticed that such large negative 
values of mayJiidxice the presence of color breaking minima at zero or 
finite temperature HLj. Demanding the absence of such dangerous minima, 
the one loop analysis leads to an upper bound on the lightest CP-even Higgs 



GeV. This bound was obtained for values of 



of 



mass of order 8C 
order (80 GeV) 2. 

The most important two loop corrections are of the form (f)^ log((/)) and, as 
said above, are induced bffi.the Standard Model weak gauge bosons as well as by 
the stop and gluon loopsEl. It was recently noticed that the coefficient of these 
terms can be efficiently obtained by the study of the three dimensional running 
mass of the scalar top and Higgs fields in the dimensionally reduced theory at 
high temperatures O. Equivalently, in a four dimensional computation o£-the 
MSSM Higgs effective potential with a heavy left-handed stop, we obtainES 



V2 



log- 



A 



H 



(304) 



327r2 



51 
16 



.9-3 



/ 42 



hi sin 1 



A] 



where the first term comes from the Standard Model gauge boson-loop contri- 
butions, while the second and third terms come from the light supersymmetric 
particle loop contributions. The scale Kh depends on the finite corrections, 
which may be obtained by the expressions given inllS. As mentioned abpve, the 
two-loop corrections are very important and, as has been shown in Ref.lI3, they 
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can make the phase transition strongly first order even for m,7 ~ OO. Con- 
cerning t he v alidity of the perturbative expansion, the /3-parameter, similarly 
to PsM, ( ^83|) , can be shown to be given by 

2 

Tn 

/3mssm ~ — 2 ' (305) 

which leads to a reliable perturbative expansion for a value of the Higgs mass 
enhanced, with respect to its Standard Model value, by a factor ~ {mt/mh)'^- 
An analogous situation occurs in the [/-direction {U = tu). The one-loop 
expression is approximately given by 

Vo{U) + Vi{U,T) = {-ml + juT^) ~ TEuU^ + (306) 

where ju and Ejj were given in El. 

Analogous to the case of the field cj), the two loop corrections to the U- 
potential are dominated by gluon and stop loops and are approximately given 

by 

log (^) (307) 

where, as in the Higgs case, the scale Ajj may onlyJae obtained after the finite 
corrections to the effective potential are computedlla. 

Once the effective potential in the (j) and U directions are computed, one 
can study the strength of the electroweak phase transition, as well as the 
presence of potential color breaking minima. At one-loop, it was observed 
that requiring the stability of the physical vacuum at zero temperature was 
enough to assure the absolute stability of the potential at finite temperature. 
As has been first noticed in ReLO, once two loop corrections are included, the 
situation is more complicatedllS. 

At zero temperature the minimization of the effective potential for the 
fields (f> and U shows that the true minima are located for vanishing values of 
one of the two fields. The two set of minima are connected through a family 
of saddle points for which both fields acquire non- vanishing values. Due to the 
nature of the high temperature corrections, we do not expect a modification 
of this conclusion at finite temperature. 

Two parameters control the presence of color breaking minima: m^, de- 
fined as the smallest value of mu for which a color breaking minimum deeper 
than the electroweak breaking minimum is present at T = 0, and , the 
critical temperature for the transition into a color breaking minimum in the 



V2{U,T) = 



167r2 



100 



^5s-2/i?sin2/3(l-4- 
9 \ mi. , 
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[/-direction. The value of m'fj may be obtained by analysing the effective 
potential for the field U at zero temperature, and it is approximately given 



Defining the critical temperature as that at which the potential at the 
symmetry preserving and broken minima are degenerate, four situations can 
happen in the comparison of the critical temperatures along the 4> (^c) and 
U (Tf ) transitions: a) < T^; fhu < m^; b) < T^; fhu > ml,] c) 

> Tc] fhu < mfj] d) T^^ > T^, fhu > fhlj. 

In case a), as the universe cools down, a phase transition into a color pre- 
serving minimum occurs, which remains stable until T — Q. This situation, of 
absolute stability of the physical vacuum, is the most conservative requirement 
to obtain electroweak baryogenesis. In case b), at T = the color breaking 
minimum is deeper than the physical one implying that the color preserving 
minimum becomes unstable for finite values of the temperature, with T < Tc- 
A physically acceptable situation may only occur if the lifetime of the physical 
vacuum is smaller than the age of the universe. We shall denote this situation 
as "metastability" . In case c), as the universe cools down, a color breaking 
minimum develops which, however, becomes metastable as the temperature 
approaches zero. A physically acceptable situation can only take place if a 
two step phase transition occurs, that is if the color breaking minimum has 
a lifetime lower than the age of the universe at some temperature T < Tc O. 
Finally, in case d) the color breaking minimum is absolutely stable and hence, 
the situation becomes physically unacceptable. 

Figure ^ shows the region of parameter space consistent with a sufficiently 
strong phase transition for At — 0, 200, 300 GeV. For low values of the mixing. 
At ^ 200 GeV, case a) or c) may occur but, contrary to what happens at 
one- loop, case b) is not realized. p-For the case of no mixing, this result is 
in agreement with the analysis ofO. The region of absolute stability of the 
physical vacuum for At ~ is bounded to values of the Higgs mass of order 95 
GeV. There is a small region at the right of the solid line, in which a two-step 
phase transition may take place, for values of the parameters which would lead 
to w/T < 1 for T = Tc, but may evolve to larger values at some T < Tc at 
which the second of the two step phase transition into the physical vacuum 
takes place. This region disappears for larger values of the stop mixing mass 
parameter. For values of the mixing parameter At between 200 GeV and 300 
GeV, both situations, cases b) and c) may occur, depending on the value of 
tan/3. For large values of the stop mixing. At > 300 GeV, a two-step phase 
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115 125 135 145 155 165 



mr (GeV) 

Figure 9: Values of ruh, rrif for which v{Tc)/Tc = 1 (solid line), = Tc (dashed line), 

mu = (short-dashed line), for mg = 1 TeV and At = 0, 200, 300 GeV. The region on 
the left of the solid line is consistent with a strongly first order phase transition. A two step 
phase transition may occur in the regions on the left of the dashed line, while on the left of 
the short-dashed line, the physical vacuum at T = becomes metastable. The region on the 
left of both the dashed and short-dashed lines leads to a stable color breaking vacuum state 
at zero temperature and is hence physically unacceptable. 



transition docs not take place. 

All together, and even demanding absolute stability of the physical vac- 
uum, electroweak baryogenesis seems to work for a wide region of Higgs and 
stop mass values. Higgs masses between the present experimental limit, of 
about 90 GeV, and around 105 GeV are consistent with this scenario. Simi- 
larly, the running stop mass may vary from values of order 165 GoV (of the 
same order as the top quark mass one) and 100 GcV. Observe that, due to the 
influence of the D-terms, values of to~~ 165 GeV, At ~ and mh ~ 75 GeV, 
are achieved for small positive values of mu- Also observe that for lower values 
of ruQ the phase transition may become more strongly first order and slightly 
larger values of the stop masses may be obtained. Remind that these results are 
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based on thc_jtwo-loop improved effective potential. Recent non-perturbative 
calculations O confirm the validity of these perturbative results. 
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